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1. INTRODUCTION 

In an effort to convert wind power to electricity in the 
most efficient way, large, slender, high aspect ratio blades 
are being used. Aeroelasticians for the last twenty years 
increasingly have developed more accurate and complex ways to 
model the aeroelastic response of such wind turbine blades. 

The overall structural benavior of such turbine blades 
was generally accepted to be sufficiently represented by an 
Euler-Bernoulli beam. Since good design often requires non- 
uniform properties and complex, twisted geometry, the overall 
equations of motion of a wind turbine blade cannot be obtained 
easily from first principles. A systematic way that this has 
been done in the past is by variational methods; i.e., Hamilton's 
Principle and the Principle of Minimum Total Potential Energy. 

It has to be noted that research done for wind turbine blades 
parallels greatly research done for helicopter blades, since 
the real difference between the two structures is only the 
structural stiffness. (Wind turbine blades are stiff er than 
helicopter blades.) 

Houbolt and Brooks in 1956 (Ref. 1) presented for the first 
time the full linear differential equations of motion for rotat- 
ing nonuniform helicopter rotor blades. They also presented a 
total potential energy functional whose variation would provide 
the same differential equations of motion. Finally, they pro- 
oosed _■* solution based on modal methods (Galerkin and Lord 
Rayleigh's method). 
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Hodges and Ormiston in 1976 (Ref. 2) presented the stability 
ana ysis of uniform untwisted cantilever rotor blades for hover- 
ing flight. The general nonlinear equations of motion were linear 
ized about the equilibrium operating position using Galerkin's 
method . 

Kottapolli and Friedmann in 1979 (Ref. 3) showed the non- 
linear differential equations of motion with periodic coefficients 
fc • a horizontal axis wind turbine. Later Friedmann (Ref. 4) and 
Straub a*.d Friedmann (Ref. 5) solved these nonlinear differential 
equations of motion by a local Galerkin method resulting in a 
finite element formulation. This was done for a helicopter in 
hover. 

The stability of a rotor blade in hover was also examined 
by Stephens and Hodges (Ref. 6) in the same year, 1980. They used 
a mixed deflection and force, formulation and they solved the non- 
linear static equilibrium equations by a collocation method. 

Kata (Ref. 7) presented the nonlinear equations of motion 
for nonuniform, twisted, horizontal axis, wind turbine blades 
using Hamilton's Principle. 

Finally, Sivaneri and Chopra (Ref. 8) in 1981 presented the 
aeroelastic stability of a helicopter rotor blade in hover by 
finite element method based on' Hamilton ' s Principle. The flutter 
analysis was done with normal rotating modes obtained from the 
finite element method. 

In the present analysis, the linear differential equations 
of motion of a horizontal axis wind turbine blade are solved 
using the finite element method coupled with the total potential 
energy formulation of the structure. 
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According to the finite element method , the structure is 
discretized into a number of finite elements. The field 
variables within each element are interpolated with respect 
to the nodal values of some generalized coordinates, using 
piecewise continuous functions. The assembly of all the 
elements will then represent the complete structure (Ref. 9). 

With the finite element method, complex geometries and 
nonuniform properties can be easily represented by averaging 
them inside each element and then decreasing the size of the 
element. 

The matrices involved in the global equations of motion 
of the assembled structure will have smaller bandwidth as 
simpler finite elements are used (e.g., the 2 node beam element 
used in this analysis) . 

For a cantilever blade, a small amount of finite elements 
is sufficient for a good displacement solution, while a larger 
number of elements is required for a good stress solution. 

In Chapter 2, the static and the vibration analysis of a 
straight elastic axis (rotating) cantilever blade was done 
using finite element method. The flutter analysis was done 
using a modal method, the modes having been the normal rotating 
modes of the blade obtained by the finite element method. 

In Chapter 3 the vibration analysis of a nonstraight 
elastic axis beam was performed. The blade was discretized 
into a number of straight elastic axis finite elements suf- 
ficient to describe the curved shape of the beam. The stiff- 
ness and mass matrices were first obtained along the local axes 
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of each individual element and then they were transformed along 
a set of axes common to all the elements. The problem was 
then solved in a similar manner with the straight elastic axis 
case. 
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2. STRAIGHT ELASTIC AXIS ANALYSIS 
2.1 Mathematical Model 

A global orthogonal system of axes, fixed in space, 

X , Y , Z is defined as shown in Figure 1, with the blade rotating 
at a constant angular velocity & about the Z axis. 

The blade is set at an offset from the Z axis, as 
shown in Figure 1. A new orthogonal system of axes x,y,z is 
now defined, with the x-axis being the elastic axis of the 
undeformed blade. The blade is also inclined to the XY plane 
by a precone angle ji and a droop angle £. The cross-section 
of the blade is symmetr ical about the major principal axis 
and Dretwisted at an angle . The shear center E.A., the 

** 'O 

center of gravity C.G., and the centroid of the cross-section 
are not necessarily coincident. 

All blade properties are in general nonuniform. The 
blade is allowed to deform under the action of rotary, vi- 
bratory and aerodynamic loadings, according to linear small 
deflection engineering beam theory. Since the blade is of 
moderate to high aspect ratio, plate bending effects are not 
significant and the deformation due to shear is neglected. 

Also torsional warping has been neglected. 

The deformation mode of the structure consists of three 
translations u, v and w and one cross-section rotation ^ as 
shown in Figure 2, under the assumption that the cross-section 
shape does not change throughout the deformation of the struc- 
ture. The three translations u, v and w describe the behavior 
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p^the elastic axis only. 

The equations of equilibrium of the blade are written 
with respect to the undeformed blade axes x,y,z. 

The aerodynamic forces are calculated according to 
strip theory based on cuasi-steady, two dimensional, incom- 
pressible unsteady aerodynamics. The position of a point A 
in the blade can be defined by blade coordinates % , and 
J which follow the blade {see Figure 2) . The position of 
the same point with respect to the xyz system of axes will be 
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where £ is equal to the x-axis coordinate of point A in the 
undeformed blade. 
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The position of the previous point A, with respect to 
the fixed system of axes X,Y,Z will be 
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2.2 Minimum Total Potential Energy Formulation 
We car. define the functional Tr where 
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Then, among all the possible displacements u,v,w and <p that 
satisfy compatibility inside and on the boundary of the 
blade, the one that satisfies equilibrium everywhere in the 
structure will also minimize the functional 7]" of equation (7). 
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£tt = o (to) 

will provide the solution for the displacements u,v,w and . 

We can express the strains in equation (8) , with respect 
to displacements u,v,w and ^ as follows. 

Let cU and els, be the lengths of an infinitesimal element 
along the blade before and after deformation respectively. 

Let AB, AC and AD be three mutually perpendicular line 
elements before deformation (see Figure 3) , with 

(d *) 1 = (ab)* +• (ac ) 1 h-(ad ) 1 


and AB ' , AC and AD' be the same elements after deformation, 
hence 

(<*0 = ( AB 9 + ( a O + ( ad ') 

Then, for small deflections 


and 
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€ » “ AB 
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Let the coordinates of point A, with respect to the x,y,z 
axes to be given by equation (1) . Then differentiating it 
with respect to ^ and substituting into the following 



dLl 

To find 


could be obtained. 

ds 
df 


, we simply substitute 

u=u“=w = ^ = 0 


into the 


ols, 

dl 


equation above. Then by substituting 

clS, dLs 

and ■■ ■." into equation 


the expressions fc * - 

oL§ 

(11), dropping nonlinear and higher order terms, we get 

£ S J = u.' - r)(w- /i cos5< +vv'sin^j -jQ' [r '*' ni £ 

+ w\oig) + OiMOtfV (1 J ) 


Substituting equation (12) into the tension equation 


T= E U 


A ‘ 


we get 

a' = — + ^(wcosjj +w"sin^) - K* f‘f 

£A 


(13) 


Substituting equation (13) into equation (12) we can eliminate 
u in favor of v, w and and hence obtain 
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Similarly from Figure 3 
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Or, substituting the expressions for AB • , AC’ and 3'C' into 


equation (15) we get 
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Differentiating equation (1) with respect to n| and 
substituting into equation (16) together with the following 
equation 



we get, after deleting higher order and nonlinear terms 


»>, = - is 4 ' 


(18) 


Similarly, substituting yj by J into equations (16) and (17) 
and reseat inq the former analvsis we get 


h = 1*' 


(«) 


Substituting equations (14), (18) and (19) into equation (8) 


for the strain energy, we get 


S-t 
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■ir"w" + [&T ♦ 

- =4 (20 ) 


where EB, and EB t are given in Append \ x M. Since the above 
equation is up to second order 


of magnitude when we come to 
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form work done by the centrifugal tension, i.e., 
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( 21 ) 


we need to consider £.|j up to second order of magnitude, too. 
Hence, if we kept second order terms into equation (14) we 
would get 

£ n = a' ♦ y[(^') 1+ (w') 2 ] ♦ tY 

- *J IT COSji 4 W*SmjJ 4 + v/iosj£)j 


+ 5 
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ir^sm jrf - w*cosj£ + ^(u-"cos^ + VY"sind) 

^ '0 


(22) 


Then, substituting equation (22) into equation (21) and ignoring 
the u' term since it is one order of magnitude higher than the 
other displacements, we get ^ 

WokK oJ ^en-fri Tension — j + (^0 


+ e /l (ir / cos^ + w*sin^) - C A - \v"c.o$£) 

-K;[j (13) 
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Substituting equations (23) and (20) into equations (9) and 

T* 

(7) respectively, and ignoring the term rrr- which does 

c A 

not contribute in the variation of 7T , we get 

f R 

TT = y J { [EWjj * Elico*jg]0 y') % ♦ [ei* 

+ EljSin'^W *) 1 + Icoifi Sih^EI, ' £I,j V W 

+(&r ♦ EB,(* , )V)‘ - 2EB.CO *^vy 

- 2. E B a 1 

R o O J 

~ J( /t x ^){-j[( tf, ) 1 + (w') X ] ♦ *| V'co^ 

+ - e A ^|Vs itiji - w'coSjrf j 


& ,R .R 

- Uvvoix - / 

f (* 

- J <^w'clx - I (j i/'clx 


(M) 


Note that rotation is about the deformed 5 axis and hence the 


t 


work done would involve ^ rather than ^ . 

The above functional 7T is equivalent to the one given in 


Ref. 1. 
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(a) Blade Absolute Acceleration Vector 
It is essential to know the magnitudes of the absolute 
acceleration components along the x,y,z system of axes (un- 
deformed blade axes) to be able to formulate the inertial 
(D'Alambert) loadings exerted on the structure. 

A systematic way this can be done without losing Coriolis 
and Centrifugal effects is the following: 

Assume a point on the undeformed elastic axis to be 
displaced by u,v,w and then the cross-section to be rotated 
by <£ (Figure 2) . 

The global coordinates of any point A of the above 
cross-section will be, according to equation (3) , 


r ' 

X 


X 


COsSt CO <><f> 

Y 

> - $ < 

J 


SIM St COS 4 

z 

y—' 

2 


•O 


where, according to equation (1) , 




f « 


r N 

X 


§ + •* 

, c 

0 

a 

2 

► — , 

f 

£ S 
\ 

+ r , 
*** 

T 

l$J 


Matrices S and F are given in equations (4) and (2b) 
respectively. 

Substituting equation (1) into equation (3) and dif- 
ferentiating equation (3) twice with respect to time, we 
get the absolute accelerations with respect to the fixed 



OMOMN. fAOE » 
OF POOR QUALITY 


19 


system of axes: 


X - 

f + 

#« 

u. 

e 

#- % 
0 

u 

o« 

i+tc 

Y “ §• 

!i\ 

## 

\r 

k + 2 £ < 
/•v 

# 

v 

» 

w 

> j 

■ + 5 ■ 

IT 


’ ( Sf *• 2SF <• SF) 


0 

1 

$ 




- COS COS P *] 

+ LJt» e 0 | - sinSlt sinp 

0 


1 


( 25 ) 


Resolving these accelerations along the x,y,z system of axes 
we get the required components 


+ 


«» 

_T 

r •• «i 

X 


U 


T o 

'o' 

u 

a 3 

= s 

#♦ 

Y 

* — < 

IT 

• 

+ 25 S 

\r 





• O 



o 

- a «. 


L*J 


w 

J 



w 

* 4 


- t ** 

+ s s 


l r 
w 



/ T»' T * o 

<S$F*2SS F 


4 

+ & e. . 

-cos(pf£)cos£ 

0 

J 

sin (f> + S) cos p 

^ v 


( 26 ) 


where 




21 




(' ) D'Alambert Forces and Moments 

Since the components of the absolute acceleration vector 
have been found, we can now proceed to evaluate the inertia 
loading thus induced on the structure. 

Assume that a point on the elastic axis of the blade is 
allowed to be displaced by a set of small virtual displacements 
fa, ftr , £w and due to the action of the real inertia load- 
ings. These virtual displacements will create a virtual change 
at the x,y,z coordinates of any point A on the cross-section, 
i.e., Sx , , &Z respectively. These virtual changes will 

be obtained by taking the first variation of equation (1) and 
using the F matrix as given in equation (2a). The expressions 
thus obtained are listed below: 

ox. - Sa ~ 

- [►) Sl'n(55+fi)t^cos(^ + fi)]W' 
- *1 cos(^ +5 !() - $siVi(j£+?f)] Si/ 

+ - $ tos 02 + /0] k/ 


Sm = Sir + 


m 


Sl = Sw + " 3sinC^ + ^) 


The virtual work done then by the inertia loading during the 
above virtual change of coordinates Cx. , S 2 will be 


( 28 ) 
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tuai 


* t: r . 




where the minus sign of the acceleration components is to 
make them inertia components. 

Substitution of equations (27) for ix, , £* into the 
virtual work equation (28) will give 




where the inertia loadings will be 



fl * 

Substituting the equations (26) for 0L X , 0^ / CX 2 previously 
found, into equations (30) for the loadings, deleting all non- 
linear terms and all u terms, since u and its derivatives are 
much smaller than v,w, and their derivatives, we get the 
equations of the six loadings as presented in Appendix A. 
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(c) Blade Absolute Velocity Vector 

To be able to perform flutter analysis, we must know the 
aerodynamic forces exerted on the blade during its rotary 
and vibratory motion, which means that we have to know the 
velocity field induced around the blade. This velocity field 
will be different at different points on the cross-section 
of the blade. To avoid this difficulty we will take as a 
reference value for the velocity, the value at the aerodynamic 
center which is at the elastic axis in our case, i.e., at 

0 = * = o (31) 

* 

Similarly, as before, allowing a point A on the elastic 
axis to be displaced by u,v,w and then the cross-section to 
be rotated by <j) , the global coordinates of this point will 
be, using equations (1), (3) and (31), 


r \ 

X 

Y 

■ - S' 

r « 

W 

IT 

* + Co < 

COS Si £ CO$P 
sinSwt co$p 

(32) 

UJ 


w 

w J 


. Sinp . 



Differentiating the above expression with respect to time we 
can obtain the absolute blade velocity with respect to the 
global system of axes 


0 * 
X 


t* . * 

* 

U. 

0 

f 1 


-SJ/SinSit CoS? 

T 

' — S * 

« 

\f 

+ s 

r* 

IT 

, + e. , 

Sltosflt cos? 

X 


o 

w 


W 

* * 


0 


( 33 ) 
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where S was given in equation (4). 

This absolute blade velocity is due to the vibratory 
and rotary motion of the blade and it is taken positive along 
the direction of the unit vectors I J K. Hence, the equiva- 
lent air velocity to this blade motion will be equal in 
magnitude with the blade velocity but taken positive along 

-4 

the direction opposite to I J K (see Figure 4a) . 

Allowing the air free stream to have velocity components 
, and w w with respect to the X Y 1 system of axes 
and the induced velocity IT to be in a direction opposite 

to WL (see Figure 4b for positive directions of all air 

w 

velocities) , then 

Totcxe air vitociiy — B^otcLe yrtLp cilj t Free c xir vtlocii ij 




- 

luclu.ee cl ire^ocitu 

f * 1 


r \ 

a 

r i 

X 



0 

0 

Y 

# 

* + < 

v w - ' 

0 

2 

V. * 


wj 





Resolving the total air velocity of the above expression 
along the xyz system of axes, according to equation (6b) we 



T 

X ♦ U, 


■ = S 

Y + V w 

W 


. H + W w - ui 



get 
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Further on we will have to resolve these components along the 
deformed blade since we are interested in the aerodynamic loads 
after deformation of the blade. Air velocities , U*j and U 5 are 

going to be taken positive along the negative L-f L»j Lj unit, 

vectors. Hence according to equation ( 6 a) 


( 36 ) 


L « 


'rcxcliai cow^pone-oi of udocrLy 
■Lanier tiat compohc.nl oi vtlotily 

per cu£cxr component of utlociiij 


Substituting equations (33) and (35) into equation (36) we get 



where F is given in equation ( 6 a) 

/v 

Since the radial component of velocity is much smaller 
than the perpendicular and tangential ones we can neglect its 
contribution to the aerodynamic forces. Also the U. displace- 
ment and all its derivatives are much smaller than the corres- 
ponding IT and W ones. Hence all the U. terms and their derive- 
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tives can be neglected in the liyj and Uj equations. 

Presenting U»j and expliciteiv (after nonlinear terms 

have been deleted) we have 

IL z: c o%4 v + siW vv +• ^ + sinc^ IT ~ S2 sm(p>+S) • 

I '9 I o , To *\ 

cos^ vV - | ^cos^p*£^ +■ C 9 tospj fi + 

+ S)cos^ + S2e 0 co$pcos^ + V w -cos^icos^p+S)- 

(wsinj^ + ircos£) - SinS2-t(cos^ - - cosS^t • 

si^(pt 5)(sm^ + ^cos^)' + V w < -*ih£ili cos(p*-k)(vv'- 

Sin^ + w'tosfy + cosffii^cosj^ - - smffii sm (M)- 

(sin^ + ^cos$j + (W w -Ui) -sin(^S)(wsm^ +■ 

+ ir cos + c<y(>tS)(si tI £ +^coi^| 

U$ = " sm ^ ^ + cos $ Vv + ffism(P-t-S)cos^ tr + Sisin(ptS)- 

in^vv - cos ( ptS)+e,, cospj f - b2f cos (£ 

- Qc.cospsing + t r w cos2teos(P+S}- 
'('vvcos^ -• ir'sm^) + smQt(si‘n^ -J-^cos^) - cosffii- 
•Si»^p+S)^cos^ - -h \^|stViS2i cos(p + £)( - v/- 

■C°^ 1- f '$**}£) “ CosS2t(s \Y\ft +■ - sinffit • 

■Si n(fi+S)(eosg + (W w - uc)|sm(P+S) • 

(“WCos^ irsm *) t cos 


sm 

< 

+< 


Also since the blade is allowed to pitch to the free stream, 

we can define this pitch to be 6 taken positive in the clock- 

0 

wise direction. Then the rate of pitch 6 , which is along the 

IT* direction, to first order terms, will be obtained' from 
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h "slpi] 


■U 


Hence, to first order terms 


n 

+[o -w u-']fd 


k 

r 

fc) 

>]* 

[0-wCrjF + [ 

)+ 


and 

Moments 


I 

J 

KJ 


4 1 

J 


The aerodynamic loads are formulated according to strip 
theory; i.e., assuming that only the velocity components per- 
pendicular to the § axis (deformed elastic axis) influence 
the aerodynamic loads. 

The convention used for positive air velocities Ukj , Uj 
angle of incidence C( , lift L , drag D , and pitching moment 
A\ . is displayed in figure 5, according to which all velocities 
and forces are taken positive along the opposite , Lj direc- 
tions while di and are taken anticlockwise as positive 
(nose-up) . 

Assuming that the blade is vertically displaced with 

C 

velocity ^ , being positive downwards, and is pitched up at 
an angle 6 to the free stream, the total lift and aerodynamic 
moment will be respectively 

L = Lc + Lnc 

M - M'. + AW 


in terms of their Circulatory and tJon-Circulatory components. 
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Greenberg theory provides the values for Lt , L nc , Me 

'arid M for two dimensional oscillating airfoil in a pulsat- 
vc 

ing incompressible flow V , the airfoil being pivoted about 
the quarter-chord point which is the aerodynamic centre and the 
elastic axis (See Ref. (71)). 

U = -jf ctc V (K + Ve t \l) 


Uc= + fa + 1 ? ) 


Q 

€ 


Me = -if^C^yv 

m nc = - t ; L c - 

Since the blade is moving with speed In downwards then 
the equivalent air speed will be upwards where 


(40 


V, = b 


Allowing the free air to approach the blade with velocity 
V (see figure 6), then projecting V and Vj, along the KJ and J 
axes and applying the conventions for positive velocities 
stated before we get 

= V COS 6 - V h Sine 
Uj = - V h cose - Vsin« 

For small angles £ , CoS£ = 1 , Sln£^i€ and since V h < V , 



Wj s - h* - Vt 


one has 
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Also, since 


then 


U, Uj « u n 

u = \/uJ + » U*| « V 


substituting equation (42) into the lift equation (41) we get 

Lc = - -j p«-<- u(-U 5 + | e) ■ 

Uc = - j I"* 6 f(- Uj ♦ | <Q 


(43) 


The drag force D will be 


D = ipac ^U : 

2 J <* 


(44) 


Constructing then figure 7, we can write force equilibrium 
equations for forces and ^ exerted on the blade and 

being positive along the positive and directions. 

= L c Sino< - DcoSe( 


F 5 = “L c cosc< - Dsm« - L 




and 


COSo( =. 


NC 


u 

l) 


Us 

smor = U 


Hence substituting equations (43) and (44) into equations 
(45) we get 

Hi = ir‘t u !‘-T iu s - %!■«.“} 

• /cvr.. • 




(46«.) 

(4U) 
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^/The convention used in this report for positive forces 
j> ^ and moments ^ ^ along the xyz system of axes, 

is displayed in figure (8) . 

Transforming equations (46 a,b,c) from the § , ►) , % axes 
to the xyz axes using equation (6a) we get for the forces 


ft] 

. = F. 

r 

0 



UJ 




or for the complete loading 

\n<f> - u"cos^ i^j + wco5^ + l rsm^)F r 

\ - 

+ (cos£ - 

?s = M«<l* 

!j= -AV 

i,= Mw' 

Note the negative sign on 0 due to *"he convention for O 

U T.J 

adopted in figure 8. 

Substituting equations (38a,b) and (39) into equations 
(46a, b,c) and then into equations (47) we can express all the 
six loadings listed above, in the following general form 
(Loading) = (Coejr^ If + (Cod^VV + + (j£oef^tr 

+ (fcoe£) vv + (Coei)^^ + (Coe^vr +- (Cocf^w 

+ (Co ef)V + (CoJ)/ +■ (CoeQw ■+■ (Coef)V 
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where the 14 coefficients for each loading are presented in 
appendix B. 

It has to be noted that many of the terms in these coef- 
ficients are shown to have negligible contribution in the load- 
ing terms and hence they can be neglected. 
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2.4 Order of Magnitude Approximations 
Throughout the analysis the following ordering scheme 
was adopted 


Quantity 


Order of Magnitude 


e 

0(1) 

e a 

0(1) 


0(1) 

c 

0(1) 

sin ( p+S) 

0(1) 

* 

0(0) 

C</«. 

0(1) 

u. 

0(2) 

IT 

0(1) 

W 

0(1) 

K 

0(1) 

< 

.0(2) 

kL 

0(2) 

k: 

Ka 

0(2) 

0(2) 

erything else . 

0(0) 


It can be seen from the above table that 1C displacement 
is one order of magnitude higher than IT and VV . This is why 
terms including U. and its derivatives have been neglected every 
where in our analysis. 

Further on, the equations of equilibrium resulting after 
the first variation of the total potential energy of the blade, 
were kept up to second order cexms 0(2) except in the torsion 
equation which was allowed to go up to third order terms 0(3) 
due to lack of first order terms 0(1). Doing this, the stiff- 
ness matrix of the structure became unsymmetric. To preserve 
symmetry, we had to allow up to third order terms in , in the 
IT and VV equations. 
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More explicitly, for all the inertia and aerodynamic 
loadings we have : 

In the ft, loading 

0(0) terms were kept together with up to 0(2) terms 
in IT and W (all their derivatives; have been neglected) 

In the V t loadings 

Up to 0(2) terms in vr and Wand up to 0(3) terms in 
were kept. 

In the lj loading 

Up to 0(3) terms were kept in LT ; Wand (p . 

Using the above approximations, many terms could be neglected 
from the inertia and aerodynamic laodings (See appendices A and B) . 

2.5 Finite Element Formulatio n 

The blade was discretised into a number N of finite elements. 
In each element, the field variables were interpolated according 
to piecewise continuous interpolation functions. To satisfy the 
compatibility and completeness requirement for monotonic conver- 
gence to the true answer as the element size decreases (see 
Reference (9)), we had to prescribe \T , IT W , w' and <p on each 
node of the structure. (See figure 9) . The corresponding inter- 
polation functions were: 

Zeroth Order Hermitean Interpolation function H° for <6 , and 
First Order Hermitean Interpolation function H 1 for IT and W . 
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Then 


r 

U“ 

»P‘ , 

H. 2 2 


% 

W 

► zz 

2 H, 0 

< 

% v 

A. 


L2 0 H.J 







1 T = [u; u;' m w'] 
f, - [ W, w,'w, w/j 

S = lh ] 



h, = [ h: h: h: h: ] 
Ho = [h; h;j 
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Substituting equations (48) ,(49) and (50) into (24) :or 
the minimum total potential energy and taking the variation 
with respect to Cf where 

f <? T <? T <f 1 

irw jb? J 

ML N 


f = 



we get 


where 


C K, ♦ K J <i = Q 

<N« /*W 


(51) 


and 


K$ : 

Kc : 

Q : 

K s = 


Sirucluro.^ SI » AAocbr/x (A bled.) 
Geomc+nc Stiffness M a^riK ( Assembled.) 
Loading or (A«e.vvi b^ed.) 

t a L t> v Lj 

\ 

b. Lt \zl 

~ ' SW X b 4 L 3 


a</ 



'Ur 


K & = 

/-W 



U 0 l,L t 
Lv ^ U 

SYA\ , 

b i La 


Q = 


< 

< 

*U N 



iUu. 


+ Ih!i Jx+Jh/oJx 

% o . 3 >t 


/R ft //T ft ft 

v : -t I \ J - 4 * | n i ' J « * J «' t *■ 

j 
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where the Li matrices are presented in appendix C and the 
ccafficients are 


l, = El.cos'g ♦ EI,s'«*si 

= EX/ +■ El* Co *f 9 

b 5 = (El a -Er^COfgsmg 

W = -EB 2 sm^(jf) 

^ = Gl +■ EB, (%') 

t ? = C„4l«^ 

t* = -e^cosg 

bi = K* 

Substituting the reduced loading equations presented in 


appendices A and 3 into equation (51) we get 

<&>*&) i- 

/*«** 

* * K s >■ Ks * K„) j = Q ( „. d] (52) 


where 



= 

Aerodynamic Mass Matrix 


= 

Structural Mass Matrix 


= 

Aerodynamic Damping Matrix 

c* 

= 

Gyroscopic Matrix 

K* 

= 

Aerodynamic Stiffness Matrix 

K s 

= 

Structural Stiffness Matrix 

K t 

- 

Geometric Stiffness Ma.rix 
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Structural "Mass" Stiffness Matrix 
Steady Load Vector 




K, 


matrices were defined on page 56. 
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where the CL K coefficients are 

a, = w\ 

0C t ■= rvi K n, 

oc, - - me 


CX^ - 

otf = 

OCfe : 

: 

a 2 - yy\ 


°4 


- M 


mecosjtf 

Sin +• ^ 

2*wS(jecos(j3 + £)^ Sinj^ 

-mS 2 1 

S2\os (P+S)[(Ki, -Kpcos(p«-S)cos2^ 

.,m(p + S)] 

bi esin^ 


- smj^ 
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CL U = «,e& l [W(P+Q + e a cos(p*$)cospJw^^. 

a u — me SJ'[S cos 1 (p*)^ CoS (p+£)cosp] cosji 
a u = mtSia cos (j> 

a is = - m&\,„(> + S)[§ cos (p+s) + e.cosP] 

a i ‘ i - -weS2 l [§tos 1 (P+£) +■ e<,cos(f$i-£)cosp] 

oc, s = -wffi’cosjj^-^^p^smj# t e?sm(p*Q- 

• cos (M) * e e 0 s i n (^j3 -i- S^) cosp] 
a i6 = 2.weS2 sin(j^+£)cosd 

a ‘* - 2**£i(j\* a - Km.) eos CP + cos 5^ Sin j$ 

a is = 2.me^ smQ}+£^ sin^ 

a K - +■ K», cos 2 ^) cos (p + ?) 

The k. [I and coefficients relate to the aerodynamic 
forces and are listed a appendix 3 with opposite sign since they 
have been obtained from the loadings directly. Hence, in order 
to be substituted in the matrices M. C„ and appearing in 

n *** fyj ft 

equation (52) their signs have to be changed since these matrices 
have been transfered from the right-hand-side to the left-hand- 
side of the equilibrium equations. 


2.6 Static Analysis 

To perform static analysis of the structure we set 

l = 0 

i = 0 

then the equilibrium equations (52) will be 

K 1 = 9 

K = K* + K s * K & + K„ 


(* 3 ) 


and 
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The solution of equation (53) will give the generalized 
displacements ^ , where 

? T = 2. m'w, w.V, 

*u w (Um. 

To obtain the shear forces Vy and and the bending moments 
and of the blade we used the same conventions for 
their positive directions that are shown in figure 8 for the 
loadings ft ' h ' % ' % , respectively, 

y vA 

from the equilibrium equations, using force summation method 
(see Reference (1)) 


The shear forces V* Va for each element can be obtained 


l^L- 


2%. 

'ix 


h 


- -h 


(5*«) 

{SH) 


. th 


Approximating the above derivatives for the i element 

21 = 1 . - Vu 

Ox l L 

_ Vzg, ~ Vk 
Ox L 


(55) 


then, the nodal values of Vy and Yj for each element will be 
obtained from 

V„., = V„ - fj,(. (St.) 

V,.. = V,. - i 4 « 


with the condition that Vv.. . = Vi,. . = o 


^ (-free tnd) 2 (A*««.ei*4) 
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It has to be noted that in equations (54a, b) we have 

- GO Aerool w no.i*iic * a) S4ruc+ itrai 


Aerodynamic 

^ ~ 0O,, e „ dWM ,, +■ OJ 


1 1 ^i^AtrojItjnzMiC Vs * 1 S Structural 

where the "aerodynamic" and "structural" parts can be obtained 

from appendices B and A, respectively, by setting all the time 

derivatives to zero and then substituting the solution for 

(which is a function of IT , W , and their derivatives) . 

To obtain the bending moments and Mi for each element 

we use the mode displacement method where (see Reference (1)) . 

M|j = (ET, + EljSin^) W*t (jEIi -EX,)sin^cosjzSir // 

- (Te„ t EBj ' Te s ftos<f> (5?< 

M 2 = (Ell - E!f)sin^tos^ w" + (EI,sin l (^ + 

-(Te fl + EB^y) C0S ^ + TeJ^i 

and S = ~k ( 58 ) 

In this case we cannot substitute directly into the moment 
equations the solution for since they contain curvature 
terms that are not represented explicitly in . Nevertheless, 

4-h ^ 

we can approximate the curvatures for the i u element, like 


w" = 



wL - vv/ • 




lx 


L 



If" — 

Kv) 


vj, - \T- 

A 


( 51 ) 


OX 


L 



(f>‘ - 

DW 


. 1. - & 



T 

Ox 


a 



and then 

substitute 

the solution for . It 

has also 

to be 

noted that 

since 

the 

equation (58) has similar 

form as 

equation 

(54a or b) 

we can 

express T in the following 

form (for 

the 
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thi 

1 element) 


wi'-h 


where 


= T7 - t « 


{ret fcwd 


(60) 


j> = Cfi + *) 

I*. V .* x iS aerodynamic i 


rodynamic structural 

and the "aerodynamic" and "structural" components can be 
obtained in an exactly similar way with j> and j) discussed 
before . 

Finally, after all the nodal values of have been found 
using equation (60) we will have to average them for each 
element; i.e. , 

T— _ T — *1.1 T 

I - - z 

and then substitute it in equations (57a,b) together with 
equations (59) and the solution for to obtain the bending 
mements and /V 

2.7 Vibration Analysis 

To perform the vibration analysis of the structure we set 

% =o 

all aerodynamic terms = 0 

Q’SftoUj = 0 


Then equation (52) will become 


M s ‘4 + C^ s+ ^ + ~ 0 


The above equation can be solved as a generalised eigenvalue 
problem (see Ref. (10)) of the following form 
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if 


where 


Ck. = M,f A 

^ ~ J Q a mol 9 = Q 

^ ^ iw r-w ' 


A = 


u). 


X 

W x 


LU 


*1 J 


»* Ht ei^enwoc^ctc. Mi<x+ru 
( m x m ) 


and 


^ ^ $ a • • • ] •* He ei<j€Mvtfc+or matrix 

( tn xm) 

for w degrees of freedom (generalised coordinates) 


(a) Uniform Cantilever Beam 

The lowest 6 vibration modes were obtained for a non- 
rotating uniform beam cantilever, discretised into 10 elements, 
using subspace iteration. 


The frequencies corresponding to these 6 modes were also 
obtained theoretically using 



where the blade properties, the coefficients OC n and t> n to- 
gether with the theoretical and numerical vibration frequencies, 
are presented in appendix D. 
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The agreement between the theoretical and numerical 
results is very good. 

Also, in the same appendix, the vibration frequencies 
for the same cantilever, rotating with iMj = 4.191 rad/s, 
are listed, for the lowest 6 modes. 


(b) NASA MOD-O 100KW Wind Turbine Blade 
The above blade was discretised into 10 elements with 
non-uniform properties which are listed in appendix E. 

The vibration mode-shapes and frequencies for the lowest 
10 modes were obtained. 

Figures 10, 11 and 12 present a comparison between the 
first flap, first lag and second flap rotating mode shapes, 
respectively (Q= 4.191 rad/s) as obtained numerically using 
subspace iteration with our analysis, versus numerical results 
presented by Lockheed California Company. 

The agreement is very good. 


( c) Investigation of the Direction of Motion - 
- Fl at Plate 

In figure 13, the direction of motion 0 M of the MOD-O 
blade during the first flap, no'n-rotating mode shape is 
presented versus the pretwist , where 




W 

i r 


It has to be noted that according to figure 13, although 
the blade is pretwisted from the root to the tip by as much as 
26°, the direction of motion doesn’t seem to be affected a lot. 
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MOD-O First Out-of-Plane (Flap) Mode Shape 

Figure 10 
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MOD-O First In-Plane (Lag) Mode Shape 


Figure 11 
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MOD-O Second Out-Of-Plane Mode Shape 


Figure 12 
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In fact, the blade appears vibrating in most of its span along 
the direction normal to the chord at the section around 30-35% 
of the span. 

Similar behaviour is exhibited in the North Wind Turbine 
(figure 14) with total pretwist 14 °, and in the McCauley Pro- 
peller (figure 15) with total pretwist 22°. 

Examining the former three blades for the direction of motion 
during the first lag non-rotating mode shape (figures 16, 17 and 
18) no such behaviour is encountered. 

In an effort to see how far we can stretch the direction 
of motion concept, we examined the vibration modes of a non- 
rotating very thin uniform plate, which has been treated as a 
cantilever beam. (See appendix F) . The plate was given a 
very high bending stiffness in the lag direct: , with a very 

small bending stiffness in the flap direction (ill-conditioned 
problem) . The first flap and first lag modes we examined for 
30°, 60® and 90° root to tip pretwist and the directic., of motion 
was plotted versus the total pretwist in each case (figures 19, 


20 , 21 , 22 ) . 


For the first flap mode shape a similar behaviour was 
encountered where the blade appeared to move throughout most of 
its span, in the direction normal to the chord of the section 
around 15-20% of the span (figures 19, 20, 21). 

No such behaviour is encountered during the first lag mode 


shape (figure 22) . 
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North Wind Turbine First Flap Mode Shape Direction of Motion 

Figure 14 
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tu = 2M Hz 

Q = o 


First Lag Mode Shape Direction of Motion 


Figure 16 





North Wind Turbine First Lag Mode Shape Direction of Motion 

Figure 17 
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McCauley Propeller First Lag Mode Shape Direction of Motion 


Figure 18 
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Flat Plate First Flap Mode Shape Direction of Motion 


Figure 19 
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Flat Plate First Flap Mode Shape Direction of Motion 

Figure 20 
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Flat Plate First Lag Mode Shape Direction of Motion 

Figure 22 
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2.8 Flutter Analysis 

To perform flutter analysis of the structure we set 




r: 0 


and hence we have 


h- Cc^ + K<£ =0 
with M = + M s 

fv r-f 

C - C a + G g- 

^ ^ ^ 

K = K* +■ K* + K* - K m 

Doing a modal analysis, we can select a number of modes P out 
of a total of W (see Vibration Analysis) and then transform 
the above eouations into this reduced modal space, through 


the followincr relation 


to obtain 


with 


% = % ^ 

Hr ?i + Ei = ° 

m = t ; m f f 

c = $; c 

s = 


Further on, setting 


we obtain 


where 


Y =. ? 


Y = ” Y 


(Gl) 


- M"K -M' 1 C 
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Equation (62) is much more handy from the computational 
point of view due to the appreciable saving in computer storage 
that results. For instance, in a typical modal analysis of a 
cantilever blade, the assembled, damping and stiffness matrices 
are 55 x 55 in size (for 10 beam elements). Picking up 5 
normal modes; i.e., first and second flap, first and second 1- 
and first ■; orsion, (p = 5) then the resulting matrix will 
be 10 x 10 in size. 

Equation (62) can be solved as a complex eigenvalue 
problem. For 


Y = t T ««i T = % 


a 


equation (62) will be 


where 


4 A. = H t 

^ -w , — . ~ 


being the complex eigenvector matrix 


and 


\ 


Ac = 


y. 


(2j,x2t>) 



being the complex eigenvector matrix 


68 ... ^ > .<■ ' 

,, «• ova \pnimku woe is 

• ; Of POOR QUALITY 

i.e., a complex eigenvalue will be 

\ = y, t <> 

u. represents the damping associated with a given mode 
while u; is the frequency of that mode. 

When is zero we have flutter, while when it becomes 
positive we have dynamic instability. 

Three numerical cases were examined, concerning cantilever 
beams., as a comparison test between results obtained from this 
analysis and results presented in the following papers, respec- 
tively : 

Ref. (3) by Kottapalli, Friedmann and Rosen 
Ref. (6) by Stephens, Hodges, Avila and Kung 
Ref. (8) by Sivaneri and Chopra 

In all cases, the vibration analysis of the beam was done 
with the beam having been discretised into 10 finite elements. 

(a) Ref. (3) by Kottapalli, Friedmann and Rosen 
In this case the response of the NASA/DOE MOD-O WIND TURBINE 
blade was examined. The blade was allowed to experience aero- 
dynamic loads due to a constant wind velocity W^along the K 
direction, while rotating at an angular velocity ^ - 4.191 rad/s. 

Note that according to our convention for wind velocity in figure 
4b, K for this case has to be negative. 

Vibration analysis was first performed and 5 natural modes 
were obtained (first and second flap, first and second lag and 


first torsion) . 
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Flutter analysis was then done, using the above five 
modes, for the following wind cases: 


Mat b 
QUALITY 


U w - o 
V* --= 0 

and W = -IS9”, “313 , “ ^ } 0 , ~ 6 2} and - 9 ^0. 5" »**/ 

W ' * 


The real part of the complex eigenvalue of the first lag 

mode i. e., u was nondimensionalised with Q and plotted 

Ww O 

against the nondimensional wind velocity ^ , where A is 

the radius of the MOD-Q blade (see figure 23) . On the seme: 
plot, results presented in Ref. (3) are also shown. 

The agreement between the two curves is good, considering 
the fact that in Reference (3) the input used was based on data 
concerning a much earlier design of the MOD-0. The blade was 
found to be stable for all the W w cases examined. The lag 
stability was found to increase with K velocity. 

(b) Ref. (6) by Stephens, Hodges, Avila and Kung 
A uniform cantilever blade was examined, with properties 
as listed in the above Ref. (6) as "Reference [1] Configuration." 

Two distinct cases were examined, the soft-in-plane and 
the stif f-in-plane case, both for zero pretwist. The properties 
associated with these cases are listed in appendix G. 

The lowest six rotati g mode shapes were obtained from a 


vibration analysis; i.e., the first, second and third flap, 
first and second lag and first torsion. 
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Flutter analysis was then performed using the above six 
modes. The complex eigenvalues obtained were non-dimensionalised 
with the rotation speed and then compared with the ones 

presented in Reference (6) . The comparison between the first 
flap, first lag and first torsion modes is presented in appendix 
G. 


It can be seen that all the eigenvalues, except the real 
part in the lag, have very good agreement for both soft and 
stiff-in-plane cases. 

The blade was found to be generally stable for both cases. 
It can be seen also that the stiff-in-plane blade is more stable 
in flap and less stable in lag than the soft-in-plane blade. 

( c ) Ref. (8) by Sivaneri and Chopra 

A uniform cantilever beam was examined, with the properties 
given by Ref. (8) and listed in appendix H. 

The blade was examined at 


Cr 


0-1 


where 


cr 

C-r = rotor "thrust coeff i‘ci'e»vL 
a* x so£i'ditw roct-io 


which is equivalent to the blade having been set at a pitch 
angle of 12.7°. 

Vibration analysis was performed and 5 normal modes were 
obtained; i.e., first and second flap, first and second lag 
and first torsion. The 5 corresponding eigenvalues were non- 
dimensionalised by S2 , the rotation speed, and then compared 



ORIGINAL PAGE IS 
OfL POOR QUALITY 


72 


j.'.A'-i , iAv'l- .• 

v T ; » 

OF POOR 

w itn results presented in Reference (8) . The comparison is good 
and is presented in appendix II. 

Flutter analysis was then performed with the previous five 
normal modes and the complex eigenvalues thus obtained were non- 
dimensionalised by ^ and then compared with the ones of Ref. 

(8) . The comparison is good and is shown in appendix H. 

The same process was repeated for 

= o.o 

CT 

or blade pitch = 0°. 

The comparison of the o^iplex eigenvalues is good. 

When the blade was set at 

ir - o.3 

cr 

or blade pitch = 29.1° 

the results obtained failed to predict the instability in lag 
which is shown in both Ref. (8) and Ref. (2). 
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3. NON-STRAIGHT ELASTIC AXIS ANALYSIS 

3.1 Mathematical Model 

A blade with a non-straight elastic axis can bo modeled as 
an assembly of a finite number of straight elastic axis elements 

A A 

with different inclinations o< and with respect to the x, y 2 
system of axes of the previous section. 

Or more explicitly, we define again a global fixed in space 
system of axes X Y2 in an identical manner with the straight 
elastic axis case. (See figure 24) . 

We define an 2 system of axes at an offset e 0 from the 
2, axis and being inclined at angles and £> from the X Y 
plane. is the precone angle and S the droop angle (See 

figure 25) . 

A local system of axes X L is then defined for every 

element (See figure 26) . The direction of these axes can be 
obtained from the y. y * system of axes if we rotate it about 

A 

the 2 -axis by an angle o( and then about the U -axis by an 
angle ^ , the angles a' and p being different for each element. 
The X L -axis is the elastic axis of this element in the unde- 
formed blade. 

The position of a point A in the blade can be defined by 
the blade coordinates ^ ^ ^ associated with the blade element 
to which this point belongs. The position of this point with 
respect to the X L ij 2 t system of axes of this element will be 
(see figure 26) 
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N 

X L 

II 

»• * 
UL 

. + F, 

0 

1 

(«) 

.V 


W L 


U j 



where 


' Ul, W L anc * ^ are e ^ ast ^ c ax i s displacements and cross- 
section rotations of the blade element in question, U. L , 1/^. and 
W u being positive along the positive Lx c directions, 

respectively, ^ along the f direction . ^ is equal to the 

X L axis coordinate of point A in the undeformed blade element, 
and 


F 

r— 


1 -w u cos(g+j£)-w'sm(jH {) 


f 

CO 5 0 H) 


W L ' 




The position of this point A with respect to the 
system of axes will be 


r 

X 


* p 

*0 

M 

+ 

r n 


• = , 


1 

UJ 




2 U 




where 



A * 

CO$c< COS p 

* 6 

•sin* CoS p 



- sinS 


cos* 


0 


A 

-sin p cos* 

A 

stnp s»n* 
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and X 0 ij are the X i coordinates of the left-hand-side 

node of the element in question in the undeformea blade, defined 
as 


M 


fc 

i 

ii 

V 

0 

UJ 

* # - 

0 


( 6 ?) 


= element length 

The position of point A with respect to the XY % system 
of axes will be 


X 


r ' 

% 


CoS Si! t CoS ^ 

Y 

‘ - 5 < 

a 

• + e c , 

SinS-t Co$^$ 

.2. 


» «• 


sm<^ 

L /o J 



where 



cosfflt COS 

-sinScji 

-ccsSt sin(p+£) 

II 

sinSt cos(p+£) 

cosffit 

-sinSt sin(pt^) 


Sin(p+£) 

0 

cos(p+S) 

0 

aid also 

be noted that 
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The rest of the blade properties are identical with the 
straight elastic axis case. 

Finally, once the displacements U. U" W and along the 
X j £ system of axes are found, then the new coordinates 
( * m ^ ) of the nodes of the blade, the new element lengths 

new A A 

( li ) their new inclinations ( c( and ( p ) ktCtAt and the new 

rife Hr 

pretwist of the blade ( <p ) will be obtained as follows 

<o NEW 
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3.2 Minimum Total Potential Energy Formulation 
Defining the functional 7T where 
77 ~ (Strain Energy) - (Work of External Loads) 



and 

Strain Energy 



C + + G ^ dxd H di 

Jt E 4 + Gy* t Gyp ^ ^5 ( 


where V* — ^*o-6xnr>e. fhe bfcccU 
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Work of External Loads 





• I ^ j( € st d 1 d $ 

/ A A 


+ 





+ 


ft- 

Q wjc/j^ 

J l Jv 


0 



The strain equations that have been obtained in the 
previous chapter (equations (12), (18), (19)) are applicable 
here since the energy equation (72) is written as an assembly 
of straight elastic axis element contributions. 

Since wp intend to preserve U. L as a separate degree of 
freedom, we will not use the tension equation to reduce it 
in favor of ir , w and (£ but we will keep it in the strain 

L L 

equations. Hence 



<JL'i - r\(yi cosjz* + W L "$injrf) - 

'05 Off 


+ wjcos 





/ 
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Strain Energy 


where 



Eli = 

EB, = 

EB, = 


r EI,co^](ui") + [EI,cos^, 

+ F-IiSm^Kw") 1 + 2 [eTi -E l, ■ 

'Vl"w u "si vKjJcozyi + [gI + EB, • 

(f) }(#') - 2 EA [ e* cos^ i r \ 

"* *' - fC^y] - 

- 2EB 1 [u‘ l 'c 04^ + W L / ''5lKl^']^'^ , |>i5 

El, + EAe* s 

EB, * EA K* 

EB, + EA K A e A 

T c f-< -UCO'+K) 1 ' 

^r 1 ' 1 ^ J 

e * [^" c0 ^ + w L "sm^ ] - e A ^ 

^ 1 - 



- W, 


Work of External Loads = 
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fit 01 + k k + v : + + A ^ 


or 


Ir = 



*0 + [^*i i,r > 1 ^ + EIj cos^j^UL*) 

4ei i <-°ty + EliSin J ^](w L ") 1 + 

1 ‘ £T,j U['V L "sin^cos^ + 

■ [GJ+ ES.Cf) 1 ]^) 2 - 2EA • 

[ e /t CoS ^ ^ ' + - K„ 2 g ^ ' u' 

' 2 EBifui'cos^ + w/si 


* 

T 



c l~ Ul ~ .? [CO + (V) ] + e„[ic"cos^ 


- e /i ^[l/Isirt^ •• Ktoij#] 

t “K«[±(*T ff]A 
'Jli^ *k K 

~ l%t d 5\ (>5<x) 
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and 


where 


7^ = J f x dx L + l^COsScosp t TtjSMrfCOSp + 

+ Ti sin pi 


(>5b) 


(V 


* ' 

’ - T ' 

jk dx ‘ 

•T, 


\\ 


1 


tvi L* 



aiirec+V#** 


3.3 Forces and Moments Applied on the Blade 
(a) Blade Absolute Acceleration Vector 
Assume a point on the undeformed elastic axis to be dis- 
placed by LL t , ir , and then the cross-section to be rotated 
t,y . 

The global coordinates X ' *Y ' 2 of a point A on the 
above cross-section will be, according to equation (68} 


X 

J*1 


cosffifc COS^I 

Y 

= 

i + < 

SmStlt cos^ / 

3 

l*J 


s ^ J 


where from equation 


r 'i 

X 
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r 

*o 

1 

► zz > 

1 ■ 

J 

V > 


2. 


+ T 


¥ 

<JL 



and from equation (63) 
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Substituting equation (63) into (65), then into (68) 
and differentiating the result twice with respect to time 
we get the absolute acceleration at point A with respect 
to the fixed system of axes 



Resolving these accelerations along the x ^ 2 system of axes 
we get the following components 


r • 

“c 


a. 

T 

f *0 ' 

X 

, = T, 

r*. Y 

u. t 

a 3 

■ = S - 
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Y 

VL 
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[*J 


W t 

W / 


♦2S T si U 


. f 
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r \ 


^ + Tri“‘ ) CtF-2s t sif 

' 2* W u / 
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-cos(ptS)cos]$ 
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(^) 




86 


> 'o 




' ^ 


r r^G£ 

<»p6»QO«4TY 


Or more explicitly, the equations for Ot K ,CX^ and Ct* are i ' V ' 4 ' 
presented in appendix I. 

(b) D'Alambert Forces and Moments 

1 ■ 1 ■“ “ ' “ i 

In an exactly similar manner with the straight elastic 
axis analysis, we assume that a point on the elastic axis 
of the undeformed blade is allowed to be displaced by a set 
of small virtual displacements , £>vr c » and 

due to the action of the real inertia loadings. These virtual 
displacements will change the X tj j* coordinates of any point 
A in that cross-section by £>i j • respectively. These 

virtual changes can be found by taking the first variation of 
the following equation, which can be obtained by substitution 
of equation (63) into equation (65): 


* 1 
X 

"J 

) “ < 

' • 
X. 

. + 1 . 

f * 

u; 

• + T F < 

r * 

0 

u. 




K 

w J 




Ignoring higher order terms and applying the following small 
angle approximation 


P00 *» Qumjty 


* A 

Sx = cos cos {5 Su c -*»"*&* - smpcos*6w t + 

t cosS cosp[ -¥) ( CoS ^ “ ^ Sm }^) + 5^ Sm J? t f £ cos jf)l ^ U *' 

+ cosi cos|5[- *| (swj^ *jZ$cos<gl) - ^ (cosjzS -^sin^j|Svv t 

+ ^COStf cos H- •jwjcos^ + iju^'sinjrf + Jv^'sin^ 

t^^'cos^l + Sin$ [ >] (sw^i + jA cos^ii ) + 3( cos }^ 

~# sm #)] " smpco$5[»)(t(K^-^sin^)- $(s»ng 
+ £cos^)]j ’ (?? e 

Sij =: sm5cosp8u u + cos5Sir u - smpsin<x£w t 

+ sinJcosp [-K|(cos^ + 0cos^)]Su;' 

+ smS cosp[-i)(sm^ + ^>c os^) - 5( tos ^ -^sinpS)]£w,.' 
t smScos Fl-1 w t cos^ + ^tC # sir^ + Jwjsin^ + ^U^'cos^] 
-cosS [ K)(sm^ + ? u co *$) * 5 ( cos ^ - ^sir^)] - 
-sin£ si*i [ vj (tosjrf - ^ (sm^ ♦ ^cos^)] j (^81 

= smj}Su w *► cosj3kvv w * + smp[ -^|(cos^ - 4- 

+ ^ (sm^ + ^cos^j] Su^' +. sinp[’*^(sm^ + jrfcosj^) 

- ^ (cosj^ - J^stn^)] 8 w; + |stn j5 [-17 W L 'cos^ + Vju^sm £ 

t^tr/cos^] 4 cosp [ >j(cos^ - J^s»n^) 

- i ^cos^)]j ^ 2c 


Mat m 

*• . SmS&g 

VTUA'JQ ^' vJ • 

Then the virtual work done by the inertia loading during 
the above virtual change of coordinates Sx , will be 

SW = Ijj ((-«x)Sx + (-*,)% + (-a^kij 


ft® jElm t*f\ V; 


[(-a,)Sx r (•<^)^ 


Substituting equations (78a, b,c) into the above expression .or 
the virtual work we will get 

SW = * p 2 ^w L + <^£w/ + <l&^} 

MM A I 


«£2 citMAt o 




where the inertia loadings will be 


Px L - If (-«*^°sScos^ -UySinScosji - aCjSinjOfi.^'l^ 


= // (-oc*Mn« -oc.cosS)/^'^]* 

U A 4 

p^ = jj (j*x sm p cos 5 -hCXjSmpsmi - a* cos/Q 


- jj \“ ot x Sin * ~oc 3 cos«xyf b 


Q = ] (-k|wcos^ t 'Jvv'sin^ + 5 u ‘^ os ^)-(" cx x* 

3^1. •jo 

-cosZcosfi - (XySinxcosp -cc 2 Sinpj + [ l l( s,h ^ + 
t ^cos^) + $(cos^ - (-oC x sim 2 + a^cosj) 

+ [ *1 C co<i ^ ~^ Sm s£) " $ (sm^ +■ ( a * $ in ^ cow 

+ a- sm£stnj - a 2 cos s)W*<m 5 * 
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A - a^smScosji - ci*|c(.$ 

?*. = jf ['l( toS ^ + $( SI "£ + ? Coi fy (-a,cosicos p 

A -a^sintfcosfS -cxjsm P)fi V 1 

Substituting the expressions for <X x CXij and O-i presented 
in appendix I into equations (79) and deleting non-linear and 
higher order terms we get the equations for the six loadings 
which can be presented in the following general form 

(Loading) = (c oef)u + .(c.Q? 

(Cat) W + (poeO ^ + C^ oe 0l >v/ 

(cJ)a + (cJ)j> * + (f 0 ^/ 

C c ° J )/ + C Ct 0/ + C Gt 0n 


where 


(Loading) - , j> , ^ , \ , ; 1, t 


Since in this report we investigated only the vibration 
analysis of the non-straight elast5,c axis beam, all the damp- 
ing and constant loading terms in the above loading equations 

were not included. Thus, only the inertia and stiffness load 
% 

ing terms are presented in appendix J. 
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3.4 Finite Element Formulation 

The blade was discretised into a number N of finite 
elements with straight elastic axis. In each element the 
field variables were interpolated according to piecewise 
continuous interpolation functions. Hence, in a similar 
manner with the straight elastic axis beam analysis, we 
prescribe on each node , V/ L , W L , ^6 or 

Uu A 


\ 

r 

Ho 


1 

H, 



H, 

< 


H. 

“ ** * 

V, 



where 


i - { u - 


u = 
i*. = 


i*. <} 

w Ll < w tl WlI j 


and H. . h, are defined as in equations (50). 

»v 

Setting 


(20) 


i = ■ tAlAlA? 


( 21 ) 


substituting equations (80) and (81) into equation (75a) 
we get 


Jr = 



^ M 


(% 1 ) 


*tL U eh«uo*.U 
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k s = 

kc = 
= 


K. = k 4 + k 6 

0*0 m * ** 

sLntcf urai zb fint<s **a.iri X t* e£me*tfc 

jeowiefri’c sft$*e« wcchn* ol He i?* element 

£ool<L uec/fo*- 

[U> U* K, U lJ 


L, L, L s L, b v L a 
^ L' t Lx 

sm v u. 


K. = 


0OO0 

r-» . r* ~ 

UO t, u 

U S UU 

U. 


is given in equation (75b) 


* l* /It 

-/ Hl T T, <« + J H.\i* L 

= -JVtcI,^ + f»V* + j B< 'll/* 1 - 

fjCxMt * j k\<u- * j 

-f<TXf# * f H- 
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where the coefficients are 

b, =■ EI.cos^ + EI 2 sm 2 ^ 
b x = El, s my + EI 2 cosy 
b 3 = ( EI 2 - EI.) cosjrf smy 
- - E B 2 cosy (y ') 

W = - EB Z smy(y') 

= g j + EE.fy;) 2 

t, = easing 

= " V 08 # 
tl - 

EA 

k - EAe A cosy 
b, t = EA e A sm^ 

k s = -EAKlf 

and the L: matrices are given in appendix C, 
Substituting the loadings presented in appendix J into 
equation (82) we get 


y ! Vj i *•' * » 


TT 


= / [ K!H + 


N doMtM.it 


where W L - mass of +k* L* e.itUi^L 

— Mass s of- *H»*. L**’ t&Unt+ii 
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and 


m 


L 


E.L, ^E<L,. E,L„ Q 

E,L 4 + K t Lv K*U QxL T , 

E,L<+X,L4 0,£ 

5YAV ^©vL, 



E-a L* j + ^L»o : E-i©L* *E lv L,. 

E„L, 

1 

II 

$ 

+ K, L> t j + K* L* «■ Li, 

si. ~ ” . . 

2., C, 

♦ K, L T » 


\ 

Hio E ^ 't H*L* 

HaC, 


SY M 


,X m L t > + X„L, 


T,g,. 


the E; , Zc , Hi , 0i , K; . Xi coefficient!, are given in 
appendix J. Since the coefficients of the mass matrix are 
much simpler than the ones of the stiffness matrix K*. > they 
are listed below 

E, = -w 

E 6 - wecos^ 

Kt = -t*l ( K‘,sir 1 ^ + K’jCOS^) 

Kg = - W ( K«, ~ Km, ) COS^ 

0 2 — mesi n^. 
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0 3 - -wecos^ 

T> = -^(K^cos 1 ^ + k*,W^) 

0-, = - m K« 

From equations (70) and 71b) we have 


X' 

T 

Lx 

Lj t 

L u 

V- J 

. = J , 

3 


Hence we can write 



where UL U" W an.5 <p are nodal displacements and rotations 
measured along the xtj2 s y stem axes. 

Rearranging the terms in equation (84) and changing the 
sign of terms we get 
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' a t ' 

1 

a N 

1 



if 

1 


= T. < 

/ 

tr 

} (85) 

K 

— 

W 


K 


w 


[tj 





where 'X* is a modified transformation matrix obtained from 
equation (84) and is presented in appendix L. 

Rearranging equation (85) to include both nodes of a 
given finite element we get 

r u = t a =L («) 

where 

1 = ^ a ' u i ^ v/\r t ir t ' w, w,'w» W x ' <$ ?ij 

It has to be noted that for this particular arrangement of 

displacements and rotations in the <£ matrix, the elements 

*0 

of are related to the elements of T« by suitable rearrange- 

ment. In appendix L a differently arranged <£ matrix is pre- 

sented which leads into a simpler transformation matrix . 

0+0 

Substituting equation (86) into equation (83) and set- 
ting the first variation of 7 1 to zero we get 

Stt = + K<^ - o (Sft) 
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where 


% 

K 

M 





'i '<*’ l<V<f,,'f»n 


»TI 




- zL I r (k.* kZk.)T. 



*U N tGtm^xit 


3.5 Vibration Analysis 

Equation (87) can be solved in an exactly similar 
manner as in 2.5 Vibration Analysis, as a generalised eigen- 
value problem of the following form 

K $ = M f A 


where 




4kfc 


£<j«nwot 


l at *ux+ 




^ ^ ^ • • . f**, tern# 4k<, ei«<HV«c4©*- ***4 rU 

L — ~ ~ 4 * 


for m deg? — s of freedom (generalised coordinates) . 


A l PAG* V; 
«****>*" 
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(a) Uni form Cantilever Curved Beam 
A curved cantilever beam with its elastic axis in the 
shape of a circular arc, was discretised into 7 finite elements 

A f 

with the following inclinations d and p . 


Element 
t* u mber 

1 

2 

3 

4 

5 

6 
7 


o’ 

10 * 
20 * 
30* 
40* 
50* 
60* 


P_ 

0* 

0* 

0* 

0* 

0* 

0 # 

0' 


The free vibration mode shapes of the lowest 5 non-rotating 
modes were obtained using subspace iteration. 

The same cantilever beam but with a straight elastic axis 
was analysed and the lowest 5 vibration modes and frequencies 
were obtained. In appendix K a comparison between these fre- 
quencies and the ones calculated theoretically, is shown. The 
agreement is very good. 

Figure 27 shows the first out-of-plane vibration mode shape 
of both straight and circular elastic axis cantilever beams. 

This mode shape consists primarily of w displacement, the UL 
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• v ' 

-98 ! . 

•’ 

i 

and U" being almost zero. It can be seen that there is almost 
no difference between the straight and the circular elastic 
axis beam mode shapes, the circular one having a slightly higher 
frequency than the straight, 


i . e . , 


\X) 


Circa. ta.r 


1.86 Hi 


W*WjU = l-SV Hi 


Figure 28 shows the first in-plane vibration mode shape 
for both beams. This mode shape consists primarily of if displace- 
ment (upper plot) . The lower plot shows the U. displacement 
associated with this mode. (VVxO) It can be seen that the cir- 
cular beam exhibits a little different mode shape than the 
straight one and it has a higher frequency, too. 

i - E -' Wc.Vc.e.r = I.JVI Hi 


= l-SV Hi 

Figure 29 shows the second out-of -plane mode shape for both 
beams. It consists primarily of w displacement (u.«irxO ). 

It has to be noted that the circular beam, although it does not 
exhibit a big change in the W mode shape, has a lower frequency 
than the straight one 

WcirtuW = ^ Hi 

U>*w,M = ll. S3 Hi 

Figure 30 shows the second in-plane vibration mode shape 
for both beams. It consists primarily of U- displacement (lower 
plot), lb 3 up^er plot shows the IT displacement associated with 
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First Out-Of-Plane Mode Shape 


Figure 27 
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this mod's (w*0). It can be seen that the circular beam has 
a considerably different tr mode shape than the straight one 
and also a lower frequency 

i ' e " '-"ci-ca^ = 1001 Hi 

— II.S*V Hi 

Finally, figures 31 and 32 present a picture of the vibration 
pattern associated with the first in-plane and second in-plane 
mode of the circular beam, respectively. 



Hi 


Undeflected Position 



Second In-Plane Mode Shape 
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In this report, a straight elastic axis blade was repre- 
sented by an Euler-Bernoulli beam, with generally non-uniform 
properties. The Total Potential Energy functional IT was formu- 
lated in terms of in-plane and out-of-plane displacements and 
cross-section rotations, based on linear theory of beams. 

The blade was allowed to rotate and vibrate in a uniform 
free air stream in general. The blade absolute acceleration 
and absolute velocity vectors were found from which the linear, 
inertia and aerodynamics loads were calculated. 

The blade was discretised in a sufficient number of finite 
beam elements along which the field variables were interpolated 
according to the Finite Element Method. Static, Vibration and 
Flutter Analysis were then performed by minimizing the functional 

TT • 

In Vibration Analysis, a uniform cantilever beam was 
analysed and the mode shapes and frequencies were compared with 
theoretical ones. The agreement was very good. The lowest 
three mode shapes and frequencies of the NASA/DOE MOD-O blade 
were obtained and compared with results from Lockheed-California 
Company. The agreement was good. The direction of motion of 
the NASA/DOE MOD-O blade, the McCauley Propeller and the North 
Wind Turbine blade were investigated for the first out-of -plane 
and first in-plane modes. All the three blades were found to 
vibrate, during the first out-of -plane mode, along most of their 
span in a direction normal to the chord at the section around 
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30-351 of their span. This behaviour was not generally 
exhibited during the first in-plane mode. Similar be- 
haviour was encountered on a set of variable-pretwist 
flat plates treated as cantilever beams and vibrating in 
the first out-of-plane mode. 

In Flutter Analysis , three comparison sets of cases 
were examined, involving References (3) , (6) and (8) . 

The agreement was generally good. 

In this report, the Vibration Analysis of a non-straight 
elastic axis blade was also formulated, based on linear theory 
of beams. The blade was modeled as an assembly of straight 
elastic axis elements with different inclination with respect 
to a common system of axes. The Vibration Analysis was then 
reduced in a form similar to the straight elastic axis case. 
Two, non-rotating, uniform cantilever beams with identical 
properties, one having a straight elastic axis and the other 
having a circular arc one. were analysed and their lowest two 
mode shapes and frequencies were compared to each other. It 
was found that there was very little difference in the mode 
shapes between the two cantilever beams for the first and 
second out-of-plane modes, but generally more difference for 
the in-plane modes. The frequencies were, in general, almost 
the same between the two cantilever beams for all the modes 


examined. 
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Appendix A - Inertia Loadings 


Second order terms are shown deleted. 


^ = 2 Sirvi - ,2 Sihy5 - 

- l{)2 msm(p+£)cos ( P + £jw - Q md Sjjr^+£jcos^P*£j- 
•cos^ji + mSi § cos(p+S) + miS2 C»cos (p-H^ ) 

- »vi2e.sm(p^Jcos(p+4js/h^ - wS2d cos*(p+S) * 

• cosjzi tr / - mS2 d co$'(p+£)sin^ w 7 + ruepefjji i/“ 


+ rue. 


t 55 


-mlr + mesm^ <f> + yy + 2 imS2 

£sm(j}+£jcos^ ^ + hh$2 it - wS/dsih^^ 

S^Ctos^ +* cos ^P+^^ ircoSj^ + ^ s,n ^j 


+ >n 


fc = 


-Mw - hi ccos^^ - 2mSUm (P*W * 

+ 2mS2dsm(P+f)sih^ ^ + 

+ mdS^jfr^'p+£)cosjrf ^ - hi&^siVi (p+£)cos(p+£)§ 
- mS2 Cc sm^p+Jjcosp - miS2csj>i^(p^j S,M 5^ +■ 
t meSliinj^p^Sfios (p + ^)[u*co$^ + Wsin^j 






, ^ 3 * 
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wesm^tr - mecos^h? - w K w ^5 - 2*M&e- 
sin^ji+S) sinjtf W - 2.mSi e Sin (p+C)cosj£ IT " 

- * K«,)co5 X (p^)cos2^ * 2W* 

•sin^ir + M^ejjMr^p+^Jcos^w + mffief • 

• iltt( p+^)cos(P+^)siVl^ ^ + M&ee^Sip^P+£) COSp- 

sin jzS - hi £2 ( K* a ■ K*, ( ) coS l (ptS)sm^CoS^ 

- m S2 C ^ si’m (H) cos (p + fjcosrf - rh fi*ee # • 

<.in(p + ^ ) COS jicos^ - wQe COs(p+f)[ § Cos(p4S) 

+ e.cosp][wcos^ - irsm^j + 2wQsm(p+£)- 

■ io^j£+%jZir\<p c os f a { K*. - «»,)*' + • 

• 5iV^^fJce>s(P+f)[K^ l ( cos*^ - 5#r,J ^) * 

+ K w ,(sm^ - cos^)]vr' - 2mS2 ( K mi - K^,) • 

• cos ( P + Sj COS sin (?■ - 2mS^ ( sm 2 ^ + 

+ Kw, cos 2 ^ ) cos (p*<?) W 

- 2. tviSJc cos (p+f)sinj^ it + 2 mffi cos (j3 +£) ' 
•(KijWjj t K*/c oty) e( fcoj^P+O 

+ £ # cos(p+S*)cosp)cos^ + mQ^*rf(f+&)' 

5(P + S)$lVi^W + 2m£j ( K Nx “ 

(p+<f) *»hj£cos^ - wiSe ( § cos 1 (P^f ) +• 

+ Co COS (P*£) cos pjsin^ + "iS^( K*n a si«y + 

+ cos 2 ^$ ) ^sitTj^^-f^^"co5 ^ p ■*" + mb^ 2 (Ki t 


cos 

cos 
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Q = - 2.wS2ecos(p + S)cos^i/ + - lC,V 

• COS (p+f ) Smj^COS^ ^ + mo& £ ( §COS 1 (p+J) 

• e 0 cos(p+S*)co5p) sm^ f6 + I *St\s»r(fi+S)- 
•cos^+^VoS^W + 

• Cos(P+^) COS ^ <f> “ wS2c.(f cos 1 (H) + 

+ e.cos(pt^) coi p)cos^ + 

•si M (p+£)cos(p+£)sm^ cosj^ + vnS2 2 ( “ 

- + mSi ( K* a * 

• cos l d + Km, sin^) co£ffit£f[r' - m ( K*. - 
~ KwJ s\j^Co$ijt \n' - m ( S| ^ x ^) 
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Appendix B - Aerodynamic Loadings 


Second order terms are shown deleted. 


Let 


4 = 


l = 


:os(ji+£)sm^ § * + *U^{sm2i’ 

- cos cos t ) 4 v "( -cos Sii find - 

- SIM cos *) * (x, - Ujj co$(p+S)tos^ 

S2cos(?+^J cos ^ ^ * Q c.jLpyfrcoijg + V w (-si»i£2i- 
cos<£ - cosffii>w , (|J+£ Jsinjg ) + Xv (cosS2t cos^ 

- SIM Qi^n^tftind.) T (W w -Uj) cos 

= Vw cosffii cos(p-t^Jsin<^ + siVi£2tcos(j5+S)s»Vi^ 

/ % * % * 






n 

■+ V„ Si m Sir os (p+S)coSjd 


* 

= V w Cos Sit COS ( p+6 ) COS ^ 4 

4 (w w -n) tm^SJcos^ 

CosSitsiMJ^ t SlH Sit sjW^Ptt) 


F. = - 


£ 

4 


Vw- £2 ( S i*iS2t «i‘m£ - cosSitsAi^^S) COSjg )| 
- 4 t £ &2^i*ffj3+&) 

Ft = 4* ~ ^ 


frt tik 1 1 ^'■■J Kt£Ni.wrt..LLt CtyAHB 


then 
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where - — ■ 

•=■ -c os<£ + s»Vi^ ‘f» 

Ai = — COJ ^ ^ 

B,y ♦ BaW ♦ V - + %W ♦ 

+ B 6 $ + B> it + B x w + 0 h- B.o^ + 

4 B„ w +■ i- B^vy' +■ ?iv 

where 

B, = - srn^ 

£>» = -^ eo *i{) 

B, =r ?/><“ (|) 2 

=• COi^g j^>ac[-2sm^^, + j 2>>f(ptS)- 

• •Jinjg - 2coig ~4j 
- SlVl^ j( >a< "'- ^ £2sj»£jJ+f)coS^ -(l 

+ /3rK‘^ s, '"£ + 4 coi £)] 

By = COS£}j*tt.[2 COS ^, - | fi>Vr(P+S) ' 

’ tr 'S* - — <J 

- sin^ r &2si*»fj5+J)siVij£ — 

- (l4^f)(L«Sj< + ^I SI '»'^)] 

B t = tosjg ij>ac(-|^,) 

"lM Jf at ( L) 
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Jj>ac[ 2S2 ?1 „(ptS)cos <f'l, - jS&'sirtf&Sjcosft 


J. ^ ol la i 

~ *' w ^ l/ ,a H ~ + £,smg) 


Bo = 


+ ! &* 



SlW 


Be = 


B, 0 

B„ 

Bn 


r J 

c0$ ^ P occ l 2C2s/n(ptS)si'n^ £ 2 ^m J '(/ 5 +X) • 

• sm^ + S2su^ff+X ) c°s^] 

~ S>n i£ 2f ac [ ’ - ^cos^) 

cos^ ^f«c[- 24 ^ + | 4 a«v.(P*rj -2&44J 

jf ac [ " 4 - sinffi^«i*i^+Sj- 

cos 5* + 

+ COsS&t Sln(p^S)sJH^) ► - (i+ ^L«J . 

~ sl '"g 5 - c0 *)£ ^ 

= Ipoc (- ^yj 

= -SiM^ i f «c[|fc + J|}%2] 

>ac 


— cos 


_J( ^ac[2^ - + ^2^,] 

~ Sll 1 j ^ 7 f ae [ * ^ S2r -V* Sin2ii C Os(jl+%) SlW + 

+ VwcosEtcos(?+S)sin^ 

- | 
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'IV 


coi i j (>«■<- 

-^ir c 


irt-Hy, + jS ~ 

+ 2 % %] 

Sl|n ^ jf“ c l ' ^ Sj^T^si'nSi-t - ^cosSt^ cos(p+£ 

•cos^ T + 4j/) - 

-f S 2 >*rf?£% + |S> 4 ] 

S - ft - %^] 

^ ^ ^ p-*v ( + si’nS2‘t • 

• si*(p*£) coscS ) + Vy* ( smS2 i Si W - 
- Oo$Si*t Sln(^+^)cOS^^| - (l + ^2^j! 

■ + y S2 sim(P+S’J^ 2 .] 


— rr Cr + + fv^^fs-vv + ij^t 

+ r >t / + i^vv ^ .i}? + r f\r + r* # + r a \r' + 

■+ iTi vv + r+ 

where the £ coefficients can be obtained by replacing 

Cos^ <—>—*>■ sm^ 


in the 


B; 


- 


Cos^ 


coefficients . 


?« = 


+ ^5 VV + ^3 £5 + ^ S' + £ + A e 

4 A ? ir + AjW + A, 51 S 1- 4 ,ci>' +- +- A.jir' 

v ^13 W* + A, V 
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A, 

A* 

Aj 

A 


= -S' I 


V 

A s 

A< 

A? 

A* 

A<i 


o 

0 


2 f ac (tf^ 

zp ac $f <-° % i 

■itfHfl U r 


A.. . - 


’ 10 


A„ = - 
A„ = 


A, 


k l Y 


j|> ott (|) l '(-S2 sya^+S) cos £ ) 

■jf at (^) , [ " "Uw ( U>5 £2i cos£ 

•CoS^ + 1 

ir(?f£ , oi 

-i.pac(^)'[-62siKv(^]^ 2 - &(-U»SmOtt- 

•co$(p*S)siVi^ + V w COsSt CoS (£+&)• 

• S '' M 5$ )] _ , 

= - ipac^| l [- •+ S 24 - Q,(V w lltnQi. 

• cos (£t£)cos^ - cos St COS (^9* 

■WS)j)] 

= ” Y f 010 ^ " V»£2(cos£ts,V^ 

h smSt^p-r(^^)cos^ - V*r£2(siVi&£ • 

. - (^sSZts^vf^+S) cos^)j 
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Appendix C 
Li Matrices 


,L 

= / h:V 

/ M fv 


e/x = 


12 64 -12 i& 

14 -</ 

SYM \ . 


4 


h, ,,t h: j 


x — 


l 


L 


H» Ho 


0 o 

1 -i 

0 o 

-i 1 

1 -i 
-i 1 


/T / 


H- !M X = io 7 . 


36 36 -36 34 

\ t 

64 -34 -4 

Ns -34 


•SYM 


\ 


Vl 


4 


Hi Ho d* = 


t 


'1 1 

-e o 

0 £ 


bws* 
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^ ” j H* Hi dx 
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2 It 
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‘ 134 
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a 
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Uniform Cantilever Beam 
Input (beam properties) 


M = 1x10 ~ 3 

inslugs/in. 

a , 

= 

3.52 

El, = 6x10* 

lb in 2 

oc 4 

s 

22.0 

EI a =12x10* 

lb in 2 


= 

61.7 

GJ * 4x10* 

lb in 2 


* 

1.57 

= 5.0 in 



s 

4.6 


= 6.5 in 
K = 718 in 

*>Ki = 6.725X10' 2 lns f u 9 s - in; 

<1 


Results Obtained 


Mode 


l St Flap 


l st Lag 


2 nd Flap 


-,nd T _ 

2 Lag 


3 rd Flap 


Theoretical 
ID (Hz) 


2.66 



16.64 


23.53 


l st Torsion 26.85 


46.70 


Numerical 
oo (Hz) 


2.66 


3.76 


16.64 


23.53 


26.85 


46.70 


0 = 4.191 rad/s 


Numerical 
w (Hz) 



3.78 


16.75 


23.62 


26.93 


46.76 
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NASA/DOE MOD-O 

« 32 in 
o 

8 «= 7* 

5 - 0 * 


X 

m 


Eli 


El z 

GJ 

( in) 

( inslugs/in) 

(lbs 

in 2 xl0 *) 

(lbs in l x!0') (lbs 

in l xl0*) 

0.0 

0.0246 


173 


173 

132 

71.8 

0.00829 


144 


232 

150 

143.6 

0.00751 


96 


179 

72 

215.4 

0.00699 


62 


143 

40 

287.2 

0.00725 


42 


120 

25 

359.0 

0.00699 


29 


97 

16 

430.8 

0.00570 


17 


72 

8.5 

502.6 

0.00466 


8 


44 

4.5 

574.4 

0.00311 


2.3 


19 

2.0 

646.2 

0.00181 


0.4 


9.5 

0.9 

718.0 

0.00168 


0 . 2 


0.6 

0.4 

e 

e A 

EB i 


EB : 

A 


(in) 

(in) (lb 

in"xl0 

l} ) (lb 

in J xl0 1 : ) (degrees) 


-8.2 

-4.9 

4.3 


-0.85 

24.5 


-8.2 

-4.9 

6.3 


-1.50 

24.5 


-8.2 

-4.9 

7.7 


-1.65 

16.0 


-7.5 

-4.5 

5.0 


-1.3b 

10 . 5 


-4.8 

-4.1 

3.95 


-0.9 5 

6.5 


-4.3 

-3.7 

2.25 


-0.70 

3.5 


-3. 8 

-3.3 

1.35 


-0.50 

2.0 


-3.3 

-2.8 

0.8 


-0.35 

0.5 


-2.8 

-2.2 

0.4 


-0.25 

-0.5 


-2.3 

-1.8 

0.25 


-0.15 

-1.2 


-1. 8 

-1.4 

0.15 


-0.10 

-1.6 







1^ .tiSiTH’iO:. .. 

4 
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k a 

K mi 


( in) 

( in) 

( in) 

15.53 

8,35 

8.35 

14.36 

9.27 

11.77 

13.19 

8.84 

12.11 

12.02 

7.64 

11.61 

10.85 

6.05 

10.23 

9.69 

5.05 

9.25 

8.52 

4.16 

8.57 

7.35 

3 . 17 

7.60 

6.18 

2.44 

7.01 

5.01 

1.46 

7.12 

3.84 

1.16 

6.35 
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Appendix F 
Uniform Flat Plate 


Input (plate properties) 




HI - 

0.777x10“'* 

inslugs/in 

* • 

0 


El, - 

0. 262x10“ 

lbs in 2 

t - 

0 


Ei* ■ 

0.236x10 7 

lbs in 1 

\ - 

0 


GJ ■ 

0.4x10“ 

lbs in J 

n 

0.866 

in 

EB, = 

0.142x10 ' 


Km, ■ 

0.029 

: n 

Eb x » 

0 


Km t * 

0.866 

in 


R = 6 in £« * 0 


Results Obtained 


Mode 

Untwisted 
a) (Hz) 

0“ - 30* 
w (Hz) 

0° - 60 * 
w (HZ) 

First Flap 

1 

90.26 

1 

90.55 

91.39 

First Torsion 

345.4 

664.7 

1187.1 

Second Flap 

565 . 7 

508.7 

407.7 

Second Torsion 

1044 . 7 



Third Flap 

1584.3 

1519.8 


Third Torsion 

1769.7 



First Lag 



1389.0 

Second Lag 

I 

i 



Third La>g 





0* - 90* 

u (Hr) 

92.75 


326.6 

1245.6 

1738.3 
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Appendix G 


Reference (6) Test Case 


Input {blade properties) 


Q 

s 

4.191 rad/s 

ft 

* 

718 in/s 

in 

= 

7.767x10 -3 inslugs/in 

El, 

s 

529.52x10* lb in 2 

H* 

Lil 

3 

971.19x10* lb in 2 (soft-in-plane ■ 

ei 2 

= 

6.053x10* lb in 2 (stiff-in-plane 

61 

= 

205.69x10* lb in 2 


= 

0.0 

K* t 

3 

17.95 in 


= 

e c = e = o 


Z 

21.98 in 

* 

z 

P = S * 0.0 

a 

3 

2tt 

C*. 

Z 

0.01 

cr 

Z 

0.1 

C 

3 

56.39 in 

f 

3 

50.8857x10"* inslugs/ in J 


* 

-1.854 in/s (scf t-in-plane c^.se) 


S 

-0.401 in/s ( stiff -in -pla'^e case) 




rJ 


j ft 
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ComDarison of Results Obtained 



PWfifMU 
°P POOR 


P«K m 

QMuiy 


Sof t-In-Plane Case 



r/a t l *va 

Mode 

This Report 

Reference (4) 

l St Lag 

* 

-0.00670 ± i 0.6924 

'-0.0011 * i 0.7014 

l St Flap 

-0 ,3203 ± i 1.0938 

-0.3245 £ i 1.0751 

l st Torsion 

-0.3613 ± i 4.966 

-0.3622 £ i 4.9875 


Stif f-In-Plane Case 



y*/a i i “V& 

Mode 

This Report 

Reference (4) 

l St Lag 

-0.00365 £ i 1.5089 

-0.0011 t i 1.5002 

l st Flap 

-0.3233 £ i 1.0842 

-0.3246 £ i 1.0741 

s t 

1 Torsion 

-0.3613 i i 4.9661 

-0.3625 £ i 4.9888 
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Appendix H 


Reference (8) Test Case 


Input (blade properties) 

R ■ 718 in 

Qi m 4.191 rad/s 

Hi ■ 7.767x10”* inslugs/in 

El, =» 525.2123x10* lb in 2 

EI a = 6.0515x10* lb in 2 

GJ = 33.5373x10* lb in 2 

K*. = 0 

K* a = 17.95 in 

= * € o ■ o 

= 26.925 in 

= 56.3916 in 

= 2.865° 

= 0° 

= 6.0 

= 0.0095 

* 53.2859x10”’ inslugs/in* 

= 0.1 
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Comparison of Results Obtained 


original page 1$ 

OF POOR QUALITY 


For £VcT ~ 0,:i 

w/n 




t/Si t 

i W /SL 

Flutter Analysis 
Modes 

This Report 

Reference (6) 

I s t Flap 

-0.325 ± i 1.099 

-0.31448 ± i 1.10 

I s fc Lag 

-0.023 ± i 1.525 

-0.03034 ± i 1.58 

S w 

1 '"Torsion 

-0.333 ± i 2.546 

-0.35209 ± i 2.38 
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Appendix I - Inertia Accelerations 


Let 


P, = 
h = 

K = 

Then 


- POs'(^tS) COSotCOS^ 4 I| *($+&! cos (f$+&) 
cot 1 Coic( + SiVi (p+S)cos(pt£)cosfi 

cos(p-tS) cos£co$p - $!*'(&+&) Si*) ft 
*»Vi cos (p4*S)cos«r smp + Si*'[p*6 Jeoip 


a 


* 


= |Q |4 u +■ Sti COS*(p1-£j SlV\«f IT + ^ p W 

4 - Q* [ (w Smp i- LT cos^) 4 Cos'(ptS) SiV >5 • (u>S^ 

- ?*■''"& ) + ft ( + ^ cw £.)] *] 4 
+ S2'[ U (" + v'jiVij 6)+ coj'(p+S’)sm?' 

* fj ( toi fi ~ P s,rt ^)] 5 

- COS V (p+^j 4 - Si'p^ " £2 C 0 cos(p*s)cosp 

4 Si (p + S)cos(p+£) 1* 2 Si ( ** CoS (P + S) Sm2 * 
CO$jSu “ O0s(p^)cos5 U* 4 COs(fi+S) vv} 

+ 2 S[ -cosfp+SJsirw cosp (- wsm^ - i/cos^) 4 

4 CO> (P + S) co«o( JZi + COS (p4-Sjsi h p 

cos fi]n + - cos (p + £)s»t *2 CO$P ( “ wcos^ + 

4 l/smfi) 4 COS (ptSj cosi jicosj^ - tos(pt£)* 


» ' ' 
t*'t. 


ORiaNAL * 

OF POOR QUAtrn 

A a 09 I a A / <*•* , , ••/ \ 

- Sm p> U>SoC W 4 [COS* COS P ^ - w “ v COS£ j 


1 

°° • J 

V. 

A 

i r 

A J* , 


- si*jB 

cos* ^5 coj^J yj 4 


4 COiicOS P • 

( **' / °* / . / \ A «Q / £ A 

- W COStf 4 V* Si*CZ> J + SlUrf^CoSO + SinpcoSoC 


»M0? W 


-Sj si'n£ cosP u. - £2 tosi. \r + Ssmpsii 

4 [ - Smo( cos^ ( - w'sirt^ - CT cosjzf) - CoSo? ( U>S^ - 

+ SiKjp ^ cos ?C )] ^ 

+ £2 [- St^S cos p(~ VVCOS^ 4 Kf\ iVl^/ J +■ COS* • 

(siVl^ 4 <fiC0Sf6 ) 4 Si* p S rn o< ( C05^ “ jzfsiVi jj ^ 

" Sj ^ - {)£ sin* CosjZ f 4 2 S2[ ( cos (p4(?)coS<? • 

• c os p - si'll (p+S) si *p) u - cos (j$+Sj si'*5 U" - 

- ( COS (p+®j Sl*j$ COS 0 ( 4 S'Vl (p4&) cos P ) w] + 

4 2SZ [(cOs(3+S) COS^COSp - SIM (B+6J Sinp j • 

* ( - wstVi^ -crcos^J + cos^P>+-^) SiW <p %mj - 
- ( u>s(p+S) smpcosS 4 St*i (“J^SJcosp) ^ co ^] vj 

+ 2 £&[ ( COs(£+£) COSrf COSp - SI* (p + £) Sinpj ■ 

•(wcos^ + l/V*^J + cos (p-rSj Si'nrf^coS^ 4 

4 ^ cos ( p4^J S/»np COSo? 4 SI* (p+sjcos pj ^ Si’iy/J ^ 

A n *• A •• A /\ •• 

4 Slh«COSpi/» 4 COS* If - SlnpSlKrfW 
■ * 

a S / ° 0/ • V 0#/ / \ A •• . y 

t [^S*** COS P \ * w St*£ ** IT COSJ3 ] - COSrf^» Sl*J5 - 

- sunp> si’hS ^?cos^|^ + si^Scosp (- wcos^ 4 

+ U" Si*^ j - CoS ^ COS^£ 4 SI* p SI*J <f> %\V\<f> J 
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a- 


= ffii }* 3 U. - S2 smS$.n(p+f)cos (p*s)w - Vtf 

+ Si [ - k + ircos^) * sirnJ stn(li+S)cos(f)fSJ- 

(cos^ -fwf 4 ) - ^ 1 

+ OG 1_ |* (-wcoscrf + u ' s,Kl ^) + Si»n (‘P+£)cos(p’t£) 

+ |Q» Sl^ (^+6 ) CO*(£ + £) Xo + & - S2 Sin z (j$+£j. 

• 2l 0 *•' t j e„ s»V\ (f£+£) cosj^ + 2&[ Sih(p+S] • 

y\ 

• smJ'Jsp Cl + (p+<T) cosS !/■ - siin(j$+£)*i*P’ 

cos<S W ) + 2 Q>[-sin(p+£) sih 5 cosp ( 

+ u-cos^J - sin (£+£) cosJ sm^ " s»n(|J+£)* 

• s(n|$ si'nS jzScos^ wj + (£+£) St'nsf cosp- 

• (^-w2os^ + \r%lv\^) - sin(f£+S)cofP ^cos^ + 

+ sin (j3+£) smp s»n« 5 + Si*\p u 4 

-v cos p vy 4- [ smp ( - WSiVij? - O' 2.05 ^ ) + 

## a ' r a < •</ ( t 

+ ^cos^cosjB wj ■+ sin p ( - wcos j& + i/si’n^ ) - 
- ^ S\v\ft cos^ J 
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Apoendix J - Inertia Loadings 


[cos*(p+S) cos^ + psm^j - 
-weS2 x 9 c.o5*(jS + £) +■ meS2 1 sm(p+S)cos(p + &) z* Q 

- me cos (p+S)cos p +> meQ |4 § 

M S2*K Mj [-oos5cos^ + siVi|3 smtf sm^ j - wefij 

- me. £2 sin* co <>J$ 5 

K M2 [ - sm? siVi(]3+S)cos(f5 + £) cos^ - J^si'n^j 
+ me B*l sm (^tS)cos ( j£ + &) X, ~ siM a (|3+£) 2 P 
+ e o Sm(p+5)cosp + [j§] 
m Wl - cos*(p+£) sihSsm^ + ucos^j 
m Km, £Z [ cosism^ + sm^sin^cos^j 
m Km, £2* [ SlMtf SIm(P*S)cOs(|J+S)siV>^ - jU^COS^] 
-cosicos p[2 : cos 1 (P+£)sim2coj^ + 

+ +- m a l e( -Z <> Los 1 (fi+&) + 2. • 

•si'n(fS+S)tos(|SffJ J e,cos(p+f)co tp + f.S)] 

- si*5 cos^[ m S2X\ ( - co? of cos ^ + Jm5 Sin^sm^ 

- m e S ^ - me S2 sin* cos p f J 

- Sm j5 [ - m iS S‘n5 sin ( p*£) cos (P+&) coS ^ 

-mSiV^siM^ Km, ♦ sm(p+S)cos 


JNTOtTIWAllY BIM* 



pi m 
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u = 


-W(£+S)e 0 t e 0 sivi(p+r)c os£ 4 - j^f)] 

- COSo^CoSp ^ hn£2 K*t ( ( “ COS*(p + J )sim5s?V»^ + 
t - Sinrtcos|$[wQ K m> (c.os*Sinfj 

+ Smj3 Si'n* CoS^)] - SlH^[w\S2 Km, ( 5m< * SlH^p+f) 

c.os(|J+S)smj* - J*, c °s^)] 


then 


±> 

lx 


oo' 

IT 


= E,U + Ej ir -i- E 3 w + E* tjt> + £ s 
+■ Ei w’ t tyU + EfU~ * E,W + £ 10 yzi 

+ Eh it' + E,iw' 


where 

E, = 

2 ~ 
b ~ 

E* = 


- COS* cos 


i( 


mcoswcos 


p) - 


SI > 10 ? cos ft ( W Si nor 


E. = 


cos ?) - Sin p ( w smft J 

- cose? cos ft ( - m sin $ ) * sin* cos^ fi+icos#' ) 

- cosS cos? (- msinJ3 cos*) - smP cosft ( - w Sm? 
Slni) - Sin P ( M COS p) 

- CoS* Cosft ( sin* sm^ - Sm ^ cos* cos^ ) 
- smS cos ft me ^-costfSiVr^ - smp sintfcos^J 

a A 

* smft nl e Cos^ cos p 

~ c.osrt'cos ft we ( - cos 5 cos p cos - sine? • 

cos ft me ( “smPcos j$ cosy* ) - £wjB me ( -s^p 

Vj 


COS' 
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E< = 


E* - 


E* = 


- CQ3of cosft w\t ( - cos* cos ft si ) - smPcosJB - 
■ Me ( - sma’CosjS Smj£ ) ~ Siftft me ( -sinft • 

' s '^) 

- COS^COS p ( KHlSZ }C ) " Jm *COS ft [ - VmQ SiM*. 

COS^) - SlKlfi ( SZ*MU* ) 

- cos ^ cos p ( vnffi CoS z (j$+S) Snnof'j - sm?Go.sfi • 
('S2hicosiJ - sin j3 ( - iS2*his/Hf/sm ( p+F ) • 

• cos ( P+£)) 

- COS* cos !( jt . ) - SlrtS cos j&(Q*m smp- 

sm*) ~ sin ft ( - M a V,) 

: - QOS* COS ft VVI e. cos 'M Sin* 

+ h cos 5?). - si n S cos ft [ >v> c S coseJ'siVijzi * 

+ si sintfcos^] ~ sm ft [ »i C&Z (sm*- 
Sin (ft + 6) COS ( ft + S) Sm<2* -l4COS^)j 
«,, ~ - cos* cos ft ( - meSify cos^j ~~ sm*cosft - 

(kicSIj swicosft cos^) ~ Sin ft ( ~ me. • 

M. cos <b ) 

1 3 ^ A i /% 

E (1 = ~ cos* cos ft ^ "* yvi£i)^ |x. Sim Si n* cos j $ • 

( w£$ij smtfcosftsmft ) — sjnft ( - mefl* • 

t'i) 


E, 

E. 


where 


U. ^ p ^ cc, } U are given in appendix I . 
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'a 


2 oo ~7 *• ~2 00 ~~7 °T ° m/ 

, a + it* it + ZI3 w + Zl* t zf 5 lt 

+ Z t w/'f Z ? M- + Z,tr + Z. 9 w + £, o y 

+ Z„ v' + 2,2 VS/' 


the coefficients are obtained if the following 

substitutions are made in the underlined terms of the 
coefficients. 


- COSc/ 


cosp> 


A A 

- £i n/ cosjj 

A 

- Srn ft 


sin* 


A. 

-ooSot 


0 



u •• 

1 1 *° 


1 1 •• 

L) *,*■ 

1 1 


= Hi u + 

H 2 \t 

+ 

H3 w +■ 


IT 



t 1 


1 1 

/ 

| I 



+ w + 

H ? u. 

•+ 

H*i r + 

n«?w + 

Hu, y* 


+ Hn If + Ha W 7 


the Hi coefficients are obtained if the following 
substitutions are n.ade in the underlined terms of 
Ej coefficients. 


- c os ^ cos 

A 


£ 


- Sli/W Co s ft 


SiwB COScx' 

A. 

sin B 


5IWCX' 


the 
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•«/ 


= 0 , u + ©J Cr + ©J w + 0 * ■ ► Os ir 

+ 06 w‘ t 0?u + Of ^ + ©T w + 0?,o / 

+ 0„ w' + 0u w' 


where 
©, = 

0, = 

0 3 = 


ne cos^ [ cos<? smp ( cm* cos t sihi smp 
(smScosp) ~ cos^ sihfi] 
wie*mj£ + me cos^ f co sx (-sin*) + 

| (cos 5)] 


■f Sinef Sir) 


v J /o L 

• (-sihftsinS ) 

14- \ « 

- cos ft ( cos ft )] 

0^ =r - vy\ ( sin 1 ; 

g + l C, cos ^) 


9 . 


r ^ * i 2 / A 

+ cosff$ cosP sm p hi Kw ? ( tmZsinff - s/np* 
•COStfCOS^) + SlH^S<h p ^ k:,( -*cos«?sm^ 

- SlhpSIhtf cos^) - cosp 

- SiMjd ^ tos £ snip w Km, (SIMS' cos ^ + Sim| 5 • 
•cosSsih^) + sinarsihp m Ktn. ( - C os* cos^ 

+ SiMp SinJ sin^ ) - cos ft Wi Ky^ ( -sm^cosft)] 
= cosj£ [c££*^i!ip wK* a ( -coiS' cos fcos^) +- 

+ sm<?S|yip hi Km i ( - *>ntf COS ft cosjrf)] 

- [ Cos£ SlUp W Km, ( COS* COS ft SW ft ) 




14 ? 


ORIGINAL PAGE 19 
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. i e-\ <o 


+ smtfSin 
+ 


8 v* K»*, (s»n * CoS p *\r\p )] + 

A / . *% A vl 


COSjrf [ “ COSp ( -M Kj| SlUpCOS^)] ~ SlV)^ [ -Cosff 
• ( w Smp Si*^ j] 

©4 - ( -sin^Sini + COS$/ CoS^-SlMp ) W K*^ ' 

• ^ -sm^ cosof cosjB j + ^ Sin ^ cos or' + Cos ^ Z> • 

• smor’Smp ) ^ Km, ( - SIH^ ^mc?CoS^) * 

+ ( - cos^ sin i - Siki^ cos c? sin f $) m Km, ( - cos ^' 
•cosicosp) + (coS^COsS - sm ^ Sin^SM^ j- 

• mK*, ( - cosj^ sin* cosp j + cos^j ( - coifi }- 

• ^ Ki, - siVi^ ( - cosp ) 

. smp(-cos^) 

0^ = sm^ [ -smtf rvieS u - cos£ weft st’n?c©sp] 

t cos^cos* sin^ mcffi ^ - sind smj3 ivictQVn* 
swjS “ co$p weS^ |^] 

0§ ~ SiV)^™e£2 [- sin* ( cos ’(ptS)s in*?] +• cos*- 

^-cosij] + cosjjS wc£2 [ cos* sin P (cos 1 M- 

•sm5) + stnoc smp (-cosa:J -c os^(~smi- 

• sw(j2>+£) co$(pt£)] 

0^ = siVijii rneS! [-smS itj + cosor ( -sin« smP)J 
•* cos^ wg.S^ 1 [ cosocsin jS a - sin<?sm p ($w* • 
•smp) + cospa] 

0, o = Sin^ »n s»ni st’n^ cos'^p + tfjsinf + 
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+ COS 


cotd jz ) +■ cos 5 (sin^ 1 * SW^ SlHtf COS^ j] 4 

cosjzi *i K* (cos^f-S)*!** (“ CoS ^) “ 

-5'^ ft) * cos5 (cosS Costs' - Sin^SmSsm^)] 

-S Mat U 4 COjiu] - SiVip -Sin^U 4 
!$■ It J r /• »* ^ if 

+ cosS U. j 4 COS^ mKI tSU* [ COS* Sin K' i,Vl / 

*7 \ A A / 

• cos*(p+£) sin£ + j^Cos^y + sm*sinp ( s<’n^ 4 

4 sinpsinicos^)] ~ SiVt j ^ w Kw, i5* [ cos* smp • 

(cos z ^ptSJ sin* (“cosj^) - Jj^®) t sih?sih g • 

• (cosicos^ - sinpsinP - sjV^[si;i|. 

• cos? u i- smB sim* u - cosB u ~ cos^. 

f *Js A + 1 ? . 1 — * 

• [ sm|3 cos# p. 4 |Z - cosp jz j 


Slrt© u 4 cosp Ut 4 smd -Sin# »* 

• iii • ii ^ ^ A ■* 

• ( - |i | cos<^ ^ + Cos* iv» k Mj be '( si**cosp cos^z> )| 
t cos^ ^ - sin? ^ j + cos<x *» K* • 

.a f ( " Sin* cosp si*n^)| 4 cos jzf |^ C05*Cosft »i ^i»4 g * 

£2 a (-(zcos^ ) + S1H X Sin B ^ K„ x £2*( sm?cosj§cos^) 

- COSp *1 K^StTf -[ZCOS^jj - S2^[ COS* SMp . 

• m K^^2 2 ( jz Sih^) 4 smS sinp w w smx 

•Cos^sm^) - cos ft 

= -COS^jZ f SiMjll || 4 Si'rl^ rvi K^ a S2 a [“ 3mX • 


(-usincO 


4 cos pi sin* cos 


B Sin j^l 
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* * >/* ' 4 


Quiury 


t COS* Sim* • 


+ cos^ n* ( t^ 4 [ - smP (- ^ cos ^ 

•CospeoS^zsJ + COS^ wi lS5* K m . ^ CO S* ( -)* COS# • 

+ SlUjSst** ( Si'hScosficos^) - Coif ( ~ |^S<Vi0 )] 

“ Sij^ inffi K m , [ cosjtsi^p ( -^cos^ ) + smc?smp • 
•smPcospcos^ - c osJ$ ( C °S^] 


I ,o -p •« *T ** T* ®? "7“ 

, u + 1 2 w t 1 5 w t J_ f <f> + Jl s \s 

^ + Ifif?' + I}LL t I^IT + T ? W + I jo ^ 


t 


tlul/'f I ti w' 


the coefficients are obtained from the 0C coefficients 
if all the products that DO NOT include underlined terms 
are set to zero and then the following substitutions are 
made in the products which include underlined terms 
COS* SM p ' — — — COS p( C«S p 

A * A 

of S liaft •" »*■ — — Sinarcosj^ 

* A 

- COsjS ^ - *- -SlM^ 




C° S ^ 


= K,u + Kjir + K 3 w + K,p + K 5 r 

+ Ks w ♦ K ? u. + K j v + K<j w + K l0 jzs 
+ K„ir' + K, t v/ 
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the i't coefficients are obtained from the Ia* coefficients if 
all the products that DO NOT include underlined terms are set 
to zero and then the following substitutions are made in the 
underlined terms 


QOS S' CoS JZ> 
sinS sm 3 


COS S' COS 

A 

SMiJcosjS 




- CO s^> 
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Appendix K 


Uniform Cantilever Beam 


InDut (beam properties) 


EI, 

e 

r* 

G2 

i 

M 

LA 

Km, 


El, ■ 0.3043x10* lbs in* 

«. - e A - EB, . EBi - o 

C. 7324x10“'’ lbs/in 
0.2318x10 s lbs in 2 

A A 

0( m 0 


I 


P 


0. 7324x10“ * inslugs/in 
0.2968x10 7 lbs 
= 0.32 in 


$. = L = 140 in 

Inner radius of the cross section - K * .5 in 
Outer radius of the cross section * ^ * . 6 in 


Theoretical data or frequencies are based on the 
following formulae for circular tube cantilever 
beams 
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Results Obtained 


Mode 

Theoretical 
u (Hz) 

Numerical 
a) (Hz) 

First bending 

1.840 

1.840 

Second bending 

11.531 

11.535 

Third bending 

32.294 

32.327 

First torsion 

164.0 

164.0 
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If instead 


we use 


and 


then 


where 


and 
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Appendix L 

T~,T> Matrices 


of 



u., u, un^V.ir.'w, w.'w, 
a, u;'w, w,'^ a, it, u;' w, w,' ^ j 


«•<.. ^ < Wu, Wu' £ ^ IT^ Ut; Wli w L ; £ 
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‘ •' Appendix M - Useful Integrals OF M 

^ «Mory 
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